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ABSTRACT 

A full spectral survey was carried out towards the Giant Molecular Cloud complex, 
Sagittarius B2 (Sgr B2), using the ISO Long Wavelength Spectrometer Fabry-Perot 
mode. This provided complete wavelength coverage in the range 47-196 /xm (6.38- 
1.53 THz) with a spectral resolution of 30-40 km s~^. This is an unique dataset 
covering wavelengths inaccessible from the ground. It is an extremely important region 
of the spectrum as it contains both the peak of the thermal emission from dust, and 
crucial spectral lines of key atomic (O i, C ii, O ill, N ii and N ill) and molecular 
species (NH3, NH2, NH, H2O, OH, H3O+, CH, CH2, C3, HF and H2D+). In total, 95 
spectral lines have been identified and 11 features with absorption depth greater than 
3(7 remain unassigned. Most of the molecular lines are seen in absorption against the 
strong continuum, whereas the atomic and ionic lines appear in emission (except for 
absorption in the Oi 63 /im and Cii 158 /im lines). Sgr B2 is located close to the 
Galactic Centre and so many of the features also show a broad absorption profile due 
to material located along the line of sight. A full description of the survey dataset is 
given with an overview of each detected species and final line lists for both assigned 
and unassigned features. 

Key words: line: identification - surveys - ISM: individual: Sagittarius B2 - ISM: 
lines and bands - ISM: molecules - infrared: ISM 



1 INTRODUCTION 

Sagittarius B2 (Sgr B2) is a well studied giant molecular 
cloud c omplex, located ~120 pc from the Galactic Cen- 
tre fe.g. iLis fc Goldsmith|[l990l ). It is unique in our galaxy, 
being one of the most massive star forming regions, and 
has an extremely rich chemistry. Many of the molecular 
species detected in the interstellar medium have only been 
observed towards Sgr B2. Furthermore, the line of sight to- 
wards Sgr B2 crosses the main galactic spiral arms lying be- 
tween the Sun and Galactic Centre. Cold clouds associated 
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with these spiral arms are see n in absorption a gainst the 
Sgr B2 continuum emission (e.g. lCreaves fc Willi ams 1994]). 
These facts make it a perfect target for systematic spectral 
surveys. 

The Sgr B2 complex consists of three main clus- 
ters of compact H 11 regions and dense molecular cores 
align ed in a north-south direction (e.g. [Goldsmith et al.] 
I1991J). These are surrounde d by a diffuse envelope (e.g. 



Hiittemeister et aLlfigg^ ). The far-infrared (FIR) emis- 



sion of Sgr B2 is most int ense close to the source Sgr B2 
M (Gol dsmith et all 1 19921 ') where the spectrum is domi- 
nated by thermal continuum from dust and has a peak near 
80 ^m. The continuum opacity at 100 /xm is high (3.8±0.4; 
ICoicoechea et alll2004l ) , which means that only the external 
layers of the cloud can be seen in the FIR. This is in con- 
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trast to longer wavelengths, where the continuum is opti- 
cally thin and emission lines from complex molecule s in the 
hot cores are observed (e.g. iNummelin et al.l[l99^ . Thus, 
the FIR spectrum of Sgr B2 contains lines due to simple 
molecules observed in absorption against the dust contin- 
uum, as well as lines of atoms and ions in emission in the 
envelope. 

Spectral surveys have previously been c arried out to- 
ward s Sgr B2 in the mm and s ub-mm region (jPriedel et al.l 
20041: iNummehn et all |2000|. 1 19981 : ISutton et al.1 Il99ll : 
Turnerlll99ll : I Cummins et al.lll986l ) but comprehensive, un- 



biased coverage at FIR wavelengths did not become possible 
until the launch of the Infrared Space Observatory (ISO). A 
program of several wide spectral surveys, including Sgr B2, 
was ca rried out with the ISO Long Wavelength Spectrometer 
(LWS; IClegg et aLlliggd ) using its Fabry-Perot (FP) mode, 
LOS. These were very time consuming and so were conducted 
towards bright sources in order to obtain a high signal-to- 
noise in a reasonable integration time. Even so, only two 
sources were covered across the entire LWS wavelength range 
with no gaps; Sgr B2 and Jupiter. In addition, two other ob- 
jects were observed with almost complete wavelength cov- 
erage; t he Kleinmann-Low nebula in the Orion Molecular 
Cloud 1 (|Lerate et al.ll200d ) and the Galactic Centre (White 
et al. in preparation). 

Here, we present the Sgr B2 survey in full - this includes 
complete unbiased coverage over the entire wavelength range 
47-196 fim. The detected lines are mostly due to rotational 
transitions between the lower energy levels of molecules, as 
well as low energy molecular vibrations and cooling lines 
of atomic species. The main identified features have already 
been published using results from this spectral survey as well 
as a series of targeted observation s using the LWS FP L04 
mode (see ICoicoechea et al]l2004l . and references therein). 
However, this is the first presentation of the entire spectrum 
(note that the FP spectrum in fig. 2 of Goicoechea et al. 
contains only the spectral lines, with gaps betw een them). 
A pre liminary report on this survey was given bv lCox et al.l 
l|l999ll . The data presented here will be available in their 
fully reduced form from the 750 Data Archiv^B as Highly 
Processed Data Products. 

In Sect. [2] we present the observing strategy and de- 
scribe the observations. In Sect. [3] we describe in detail the 
data reduction with particular emphasis on the improve- 
ments in calibration developed for the Sgr B2 dataset. In 
Sect. |4] we present all the data with a full list of assigned fea- 
tures and summarise the results for each species. In addition, 
we present a list of unidentified features. This is particularly 
important in the 158-196 /im range that will be re-observed 
at much higher spatial and spectral resolution by the HIFI 
instrument on board the Herschel Space Observatory. The 
spectral range >57 fim will also be re-observed by the Her- 
schel PACS instrument, with spectral resolution a factor of a 
few below that of our survey, although with ~8 times better 
spatial resolution and much improved sensitivity. 
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Figure 1. The ISO LWS beam footpr i nt ove rlayed on the 100 
contours from from [Goldsmith et aL I 1I1992I') . Offsets are centred 
at the nominal position of Sgr B2 (M) at a = 17''47™20.2=, S = 
-28°23'7.3" (J2000). The contours are 0.2, 0.3, 0.4 ....0.8 of the 
peak intensity (18 kjy). The location of the ISO beam quadrant 
where spectral fracturing occurs is indicated. 



2 OBSERVATIONS 

Sgr B2 was observed as part of the guaranteed time spectral 
surveys programme ISM_V. Full wavelength coverage (47- 
196 fj,m) was achieved over the entire LWS spectral range 
using the FP mode, LOS. In addition to this, extra observa- 
tions were scheduled as a solicited proposal (SGRB2_ZZ) in 
the range 167-194 ^m. These extra observations aimed to in- 
crease the signal-to-noise in the long wavelength range. The 
first observation was carried out on 1997 March 6 {ISO rev- 
olution 476) as a test of the LOS mode before it was fully 
commissioned. The remaining guaranteed time observations 
were carried out between 1997 April 3 and 1997 April 8 with 
the solicited proposal observations between 1998 February 
26 and 1998 March 14. 

To operate the LWS at high resolution, a Fabry-Perot 
(FP) interferometer was placed into the beam in front of 
the grating (which then acted as an order selector for the 
FP). To maintain high spectral resolving power across the 
full LWS wavelength range, one of two FPs was used, each 
with its etalons and spacing optimised for a particular wave- 
band. The shorter wavelength FP, termed FPS, had a nom- 
inal range 47-70 fim and the longer wav elength FP, terme d 
FPL, had a nominal range of 70-196 (|Davis et al.lll995l ). 
Radiation was detected over the whole spectral range using 
10 detectors, each with its own band-pass filter. In each LOS 
observation both the LWS FP and grating were scanned to 
cover a wide range in wavelength. The final data consist of 
a series of many FP 'mini-scans', each at a different grating 
angle. This is in contrast to the L04 mode, where the FP was 
scanned at only one or two grating angles, giving a narrow 
targeted observation. 

The spectral resolution achieved was SO — 40 km s~^ 
across the range (see Sect. IA2I for a detailed discussion). 
Each observation was carried out with a spectral sampling 
interval of a quarter of a resolution element with each mini- 
scan repeated S times for the guaranteed observations and 
4-6 times for solicited observations. Each data point had a 
detector voltage ramp length of 0.5 s. The final dataset con- 
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sists of 36 individual observations with a total of 53.6 hours 
of ISO observing time. These are detailed in Table ICll 

The LWS beam had an effective diameter of between 
66" and 86" across its wavelength range (|Grv et al.i ,2003.) 
and was pointed at coordinates a — 17''47™21.75'', 5 — 
-28°23'14.1" (J2000). This gave the beam centre an offset 
of 2f.5" from the main FIR peak of Sgr B2. This pointing 
was used to exclude the source Sgr B2 N from the beam. 
An in depth study of the LWS beam profile has shown that 
it was asymmetric with characteristics th at vary for di ffer- 
ent positions of the source in the aperture l|Llovdll2002b[ ). In 
particular if the source was located in one quadrant of the 
aperture, there are problems of 'spectral fracturing', where 
the spectral shape does not agree on di fferent detect ors - this 
is not a problem for on-axis sources (iLlovdl [2002al ) . As al- 
ready mentioned, in order to exclude radiation from Sgr B2 
N, the brightest FIR peak was observed off-axis. For this 
reason, it was important to maintain the same orientation 
of the LWS aperture for all of the survey observations. A 
spacecraft roll angle of ~ 93° was used for all observations, 
ensuring that the source did not fall in the fractured part of 
the aperture. A representation of the LWS beam footprint 
on a 100 ^m map of Sgr B2 is shown in Fig. [1] 

In addition to the L03 observations, the whole spectral 
range was observed at lower spectral resolution using the 
LWS grating mode LOl (resolving power, A/AA ~200-300). 
This was carried out with the same pointing and roll angle 
as the FP observations. Details of the LOl observation are 
given in Table [CT] 

2.1 Non- prime data 

During each observation, the instrument settings were op- 
timised for one detector only (so that a single FP spectral 
order was scanned over the maximum in grating response 
at that detector position). This detector was designated as 
'prime'. However, at the same time, the other 9 non-prime 
detectors simultaneously recorded data. These data are also 
useful if the combination of grating angle and FP gap were 
such that an adjacent FP order was transmitted to the non- 
prime detector position. As the grating response covered a 
similar wavelength range to the separation between FP or- 
ders, this occurred frequently. 

In the course of the Sgr B2 survey, all wavelengths in 
the LWS range were observed in prime data. There is then 
a huge additional dataset (9 times larger) containing non- 
prime data. Much of this was observed with high transmis- 
sion through the instrument and so is extremely useful to 
increase the signal-to-noise ratio in the prime data and to 
confirm the detection of weak lines. This is particularly use- 
ful when data with the same wavelength can be recovered 
from different detectors. 

The non-prime dataset is also useful for deriving an 
accurate flux calibration for the survey. The next section 
describes in detail how the standard FP calibration was 
extended by including non-prime data. In particular, dark 
current including stray light could be recovered where the 
combination of the FP and grating blocked all transmis- 
sion to a non-prime detector. Although non-prime data had 
been used before bot h in L03 and L04 observations (e.g. 
ICeccarelli et al.ll2002l ). they were always reduced in an ad 
hoc way. The comprehensive calibration developed for the 



Sgr B2 survey allows both prime and non-prime data to be 
combined in a consistent and reliable way across the LWS 
wavelength range. This technique has since been used fo r 
the LWS FP survey towards Orion KL l|Lerate et all 120061 ). 



3 DATA REDUCTION 

Processing of the LWS FP data was carried out using the 
LWS offline pipeline (OLP) version 8 and routines developed 
for th e LWS Interactive Analysis package (LIA; iLim et ahl 
l2002h . The basic OLP cahbration, to the 'Standard Pro- 
cessed Data' stag e, is fully described in the LWS handbook 
(|Grv et al. I I2OO3I ) . This produced data still in engineering 
units (FP gap voltage and detector photocurrent). The basic 
calibration was extended specifically for the Sgr B2 survey, 
with additional steps developed for the LIA FP process- 
ing routine, 'fp_proc'. This routine applies the conversion 
from raw engineering units to flux and allows an interactive 
optimisation of each step. The extension to the standard 
calibration is described in Appendix (determination of 
accurate dark and stray light values, and characterisation of 
the instrumental response and throughput). The improve- 
ments made to the fp_proc routine were included as part 
of the version 10 release of the LIA software. Further opti- 
misation (e.g. removal of glitches) was carried out using t he 
ISO Spectral Analysis Package (ISAP: ISturm et al.|[l998l ). 



3.1 Wavelength calibration and mini-scan shifting 

The wavelength calibration for the FP data was deter- 
mined and monitored usin g observat i ons o f several stan- 
dard sources as described in iCrv et al.l(|2003l ). We have used 
the latest wavelength calibration coefflcients from OLP 10. 
Adopting the most conservative value for the FPS wave- 
length uncertainty gives an accuracy of 6 km s~^. For FPL, 
a systematic test of the wavelength accuracy was performed 
using CO lines observed towards Orion - in this case, the 
observed re sidual velocity d ifferences were never worse than 
11 km s-i (|Grv et al.ll2003l ). 

However, note that the standard wavelength calibration 
was determined using spectral lines within the nominal oper- 
ating range of each FP and the uncertainties quoted above 
only apply within these ranges. We have used non-prime 
data from FPL for the region <70 /^m, and these data show 
a clear error in wavelength calibration that is larger than 
11 km s~^. No correction has been made for this systematic 
shift - it is discussed further in Sect. 14.51 

The wavelength calibration of the grating is also impor- 
tant in the reduction of FP data as the shape of the grating 
response function must be removed from each mini-scan. 
The precise location of the grating transmission maximum 
in wavelength is required so that the correct portion of the 
response shape can be used. However, the requirements for 
high resolution FP observations were not taken into account 
in the original design specification for the grating positional 
accuracy. In order to provide the necessary information for 
FP observations, the grating position should have been mon- 
itored roughly 50 times more accurately. The effect on data 
calibrated with the LWS OLP pipeline is that adjacent mini- 
scans do not necessarily join together. A small shift in the 
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Figure 2. Example of prime data from one observation (TDT 47600809; grey) with the fitted continuum. The LOl grating observation 
is shown below (thick line) and also the detector RSRF (dashed line). Two places where dips in the RSRF correspond with emission in 
both LOl and L03 observations are show by the dotted lines. 



grating position can correct for this but there is no inde- 
pendent means of determining its value except for direct 
inference from the FP data. To provide the best estimate 
of these shifts, we have used the LIA FP reduction routine, 
FP_PROC, to interactively process each mini-scan. This rou- 
tine allows the shift to be adjusted until adjacent mini-scans 
show the best agreement. In the reduction of the Orion KL 
spectral survey data, an addi tional step to fit t he shape of 
each mini-scan was applied l|Lerate et al.|[200^ . However, 
for Sgr B2 the interactive shifting was enough to align ad- 
jacent mini-scans and the extra step was not applied. 

The wavelength scale of each observation was corrected 
to the kinematical Local Standard of Rest (LSR) frame, to 
account for the motions of the Earth around the Sun and the 
Sun around the Galaxy. The corrections applied are detailed 
in Table [CTl 



3.2 Glitches 

Glitches in the data were caused by charged particle hits on 
the detectors. These caused a jump in the detector volt- 
age ramp, changing its slope (and possibly the slope of 
subsequent ramps). These appear as data points with high 
(or low) photocurrent. Very large spikes in photocurrent 
were rem oved automatically in the pipeline processing (see 
iGrv et a l. 2003 ) but smaller glitches and the decaying tails 
of large glitches could not be removed automatically due to 
the difficulty of finding an algorithm that could determine 
the difference between glitches and real spectral lines. This 
meant that the data had to be deglitched by hand. 

In the LOS observations each mini-scan was repeated 3- 
6 times and this usually meant that there was at least one 
repeated scan with no glitch signature. In order to distin- 
guish between glitches and real lines, the data were plotted 
with each repeated scan in a different colour using the his- 
togram plot style within ISAP. In this way it was possible to 
make a decision on each glitch by comparing with the other 
repeated scans. The glitches that occurred in the LWS FP 
data can be described by four generic shapes and this pro- 
vided a template against which to compare the data; a single 
spike (positive or negative), a sudden jump in photocurrent 
with decaying tail, or an upward spike with sudden fall and 
gradual recovery. 



3.3 LOl observation 

One LOl observation was included as part of the survey 
with the same pointing and spacecraft roll angle as the L03 
data. However, due to the strength of the source, the voltage 
ramps for the long wavelength detectors showed non-linear 
behaviour (this was not a problem for the FP observations 
due to the low transmission of the FP etalons). These non- 
line ar effects cau se the spectrum to sag at long wavelengths 
(see lLeekdlioOol ). 

In order to calibrate these data, we discarded the second 
half of each detector voltage ramp, effectively reducing the 
length from 0.5 s to 0.25 s. Further correction was then ap- 
plied using the latest version of the 'stron g source correction ' 
and the LOl post processing pipeline l|Llovd et al.l 120031 ). 
The final calibrated spectrum shows reasonable agreement 
in the overlap region between detectors, and the sagging has 
been removed. 

Even though the grating observation required these ex- 
tra measures to correctly calibrate it, the resulting spectrum 
has a lower uncertainty in the flux than the FP observations. 
This is because the intrinsic uncertainty in the calibration of 
FP absolute flux is much higher than that of the g rating due 
to th e additional throughput correction required (iGrv et al.l 
l2003h . The uncertainty i n grating flux fo r the short wave- 
length detectors is ~10% (iGrv et al.ll2003l ). and we estimate 
that it is ~20% for the strong-source-corrected long wave- 
length detectors. 

The final LOl spectrum was only used to obtain the 
continuum level needed to determine the absolute flux of 
the atomic and ionic emission lines observed in the survey 
in Sect. 1331 



3.4 Continuum features in the FP data 

In the final FP spectrum there are several jumps between 
data from different observations and between detectors when 
viewed on an absolute ffux scale. These are due to uncertain- 
ties in the multiplicative calibr ation factors applied in the 
reduction (|Swinvard et al.ll 19981 ). especially the FP through- 
put, absolute responsivity and detector RSRFs. A good way 
to bypass these uncertainties is to divide the continuum level 
into the data to obtain relative line depth. The remaining 
error is dominated by statistical noise in the data with a 
small additional uncertainty due to the dark current deter- 
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Figure 3. The full LOS FP spectrum after co-adding prime and non-prime data using FPL for all LWS detectors (below 70 /^m only 
non-prime data were used). The positions of detected features for several species are shown. The data were binned at 1/2 the instrumental 
resolution element. 



mination (this has been weU constrained: Sect. I A l| ) and con- 
tinuum fit. 

However, within each observation (~5-10 /.im wide), 
there are some features on scales matching the grating reso- 
lution (0.3-0.6 /im) which complicate the continuum division 
process on large scales. These are probably due to structure 
in the RSRF of each detector. They contain some known 
features due to spectral lines present in the observations of 
Uranus (see Sect. IA2[) but not included in the model (e.g. 
the HD line at 112 /im). These appear in the FP spectra as 
large scale features (see Fig. [2}. In addition, sudden changes 
in the RSRF can cause large features due to the delay in 
responding to fast changes in flux caused by the transient 
response of the detectors. 

A further effect that could cause structure on large 
scales is the fact that the spectrum is built up from many 
narrow mini-scans, each of which is interactively adjusted 
in the data reduction. This could mean that spurious large 
scale features are inadvertently introduced by the shifts ap- 
plied. However, in practise, for prime data where every mini- 
scan was centred near the maximum in grating transmission, 
the freedom in the shifting process is only in their slope and 
curvature and not in the absolute flux level. Therefore, drifts 
that were not present in the original data cannot be intro- 
duced purely by interactive shifts. 

In fitting the large scale continuum, we have assumed 
that all features in the FP data on scales of 0.5-1 /.tm are not 
real and were caused by multiplicative instrumental effects 
(the only possible additive effect is wavelength dependant 
stray-light and the optical design of the instrument should 
have ensured that this was negligible). The continuum fit 
was carried out by masking out strong lines matching the 
scale of the FP resolution (even wide absorption due to line 
of sight clouds can be clearly separated from the wide fea- 
tures on the scale of the grating as they have maximum 
width ~10% of the grating resolution). The masked data 
were then rebinned to a quarter of the resolution element 
width of the grating (the grating resolution element width 
is ~0.3 lira for the SW detectors and ~0.6 /xm for the LW 
detectors). The rebinned spectrum was interpolated back to 
the original wavelength scale using a spline fit. An example 



140 






120 


_Non — prime 
- FF^L 




100 






80 




/ Prline+Non — prime\ _ 
FPL \ - 


60 






AO 






20 




Prime FPS \. " 










40 60 SO 100 120 140 160 1 BO 

Wavelength (^m) 

Figure 4. The signal-to-noise ratio achieved across the survey 
range for data binned at 1/4 of the instrumental resolution el- 
ement. The recovered non-prime data using FPL below 70 /xm 
(dashed line) shows a clear improvement over the prime data us- 
ing FPS. 



for one observation is shown in Fig. [2] with the RSRF and 
LOl grating observation. This figure shows that in this case 
two of the largest features have counterparts in the RSRF 
(in the opposite sense to the data), but also that these both 
show up as spurious features in the LOl observation. 

The method of re-binning and spline interpolation was 
used in preference to fitting a polynomial because very high 
orders were necessary to cover all features and this meant 
that spurious features were introduced, particularly in the 
regions of masked lines. 



4 RESULTS 

The full spectrum is shown in Fig. |31 and is broken up 
into narrower wavelength ranges with labelled features in 

Fig. [ml 
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Figure 5. Prime data (using FPS) in the region 47-70 /^m are shown in grey with the recovered non-prime data using FPL shown in 
black. Both are binned at 1/4 of their respective resolution elements. Note the gap in the non-prime data (64.2—64.6 ^m). 



4.1 Signal-to-noise achieved 

The final signal-to-noise (S /N) ratio achieved in line to con- 
tinuum ratio data depends both on the transmission of the 
instrument and on the strength of the underlying contin- 
uum of Sgr B2. Figure |4] shows the variation in S/N with 
wavelength based on the RMS noise in the final spectrum 
presented in Fig. IBll - ie. with a bin width of 1/4 of the 
instrumental resolution element. We have investigated the 
absorption depth/emission of the detected features relative 
to this RMS, rather than increasing the bandwidth of each 
bin to enable a detailed search for weak features. This is 
due to the other uncertainties in the data and the potential 
(partially resolved) velocity structure in the line shapes con- 
fusing any analysis of total flux using an equivalent width 
approach. 

The maximum S/N ratio achieved is in the range 90- 
125. This falls off at longer wavelengths due to the decrease 
in continuum level. At short wavelengths in the range cov- 
ered by FPS (<70 /im), the instrument throughput is low, 
as well as there being a sharp fall off in the continuum. The 
S /N in the prime data was improved significantly by includ- 
ing the non-prime observations. This is particularly useful 
in the wavelength region <70 pim where the throughput of 
FPL is higher than FPS. We have recovered and calibrated 
all useful non-prime data from FPL over the entire range of 
the survey. The resulting spectrum below 70 /xm is compared 
with the prime data from FPS in Fig. [S] and the improve- 
ment in S/N shown in Fig. |4l The non-prime spectrum was 
stitched together from many fragments observed using FPL 
with detectors SW2, SW3 and SW4. There is only one gap in 
the non-prime coverage between 64.2 and 64.6 /xm. The only 
disadvantage in these data is a slight reduction in spectral 
resolution and some systematic error in wavelength calibra- 
tion (by up to ~25 km s~^ - see Sect. 14.5) . Table [l]compares 
the S/N and spectral resolution for several lines in this re- 
gion. For all the results presented in the following sections, 
only the higher S/N data using FPL are considered. 

For wavelengths above 70 /xm, most of the range was 
covered in at least one additional observation using a non- 
prime detector and these data have been co-added with the 
prime data in Fig. IBll The non-prime data are particularly 
useful because they extend the overlap in spectral coverage 
by adjacent detectors, allowing features in the overlap re- 
gion to be independently checked. All the data presented 
here were observed using the LOS mode and additional nar- 



Table 1. The gain in signal-to-noise between prime and non- 
prime data at selected wavelengths below 70 /xm. This is shown 
as the ratio of signal-to-noise using FPL to that using FPS. The 
resolution element for FPS and FPL is also shown. 



Wavelength 


S/N Gain 
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At) FPL 


(^tm) 




(km s~^) 


(km s"'^) 


51.8 ([Om]) 


8.5 


45 


62 


53.3 (OH) 


8.8 


45 


61 


57.3 ([Nm]) 


4.8 


44 


57 


63.2 ([Oi]) 


5.0 


45 


52 



row L04 scans were not included. The L04 scans were gener- 
ally observed with a higher number of repeated mini-scans 
per point and so have higher si gnal-to-noise (although only 
around the targeted lines - see iGoicoechea et al ] |2004l . for 
details of the L04 observations of Sgr B2) . 

4.2 Assigned features 

The majority of lines that appear in the spectral range of the 
survey are due to rotational transitions between the lowest 
energy levels of simple hydride molecules (e.g. H2O, OH). 
These are generally seen in absorption against the back- 
ground continuum emission, due to the envelope of Sgr B2 
and in the case of the lowest excitation lines (normally the 
ground state transition) , due to clouds along the line of sight. 
We also detect rotation-inversion transitions of the pyrami- 
dal molecules NH3 and IfaO'''. In addition, ro-vibrational 
lines from the lowest energy bending modes of non-polar 
carbon chains occur in the survey range. In particular, sev- 
eral ro-vibrational transitions of C3 in its V2 bending mode 
were observed. Table [2] summarises the detected molecular 
lines and gives an estimate of the column density of each 
species. There are also emission lines from forbidden spin- 
transitions of atoms and ions (C 11, N 11, N ill, O I, O ill- see 
Sect. l4.4"|l . Finally, several isotopic lines of OH and H2O are 
detected - see Sect. 14.61 These are important for modelling 
the physical conditions in the source as they have much lower 
opacity than the lines of t he main isotopologues. 

As already noted bv ICeccarelli et al] (|2002l ). no rota- 
tional transitions of CO were detected in the survey (the 
lowest energy transition in the range is J =14-13 at 186 /xm, 
with upper energy of 581 K). In addition. IGoicoechea et al] 
^2004) did not detect CO anywhere in the extended region 
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Table 2. Summary of assigned lines in the survey, giving the number of detected lines, an indication of whether absorption from the 
line of sight is observed in any of the lines, and the fitted column density at the velocity of Sgr B2. The method used to calculate these 
values (fitting the ground state line, assuming local thermodynamic equilibrium (LTE), using a large velocity gradient model (LVG) or 
using a non-local radiative transfer model) and the reference is shown in the last column. 



Sp6cies 


Features 


Line of sight 


Column density 


Fitting niGthod 




detected 


absorption 


in Sgr B2 component 


& ref. 




in survey 




(cm~^) 




NH3 


21 


no 


3x10^** 


LVG [2] 


NH2 


5 


no 


(1.5-3) xlQis 


LTE [1] 


NH 


3 


no 


4x10" 


LTE [1] 


H2O 


13 


yes 


(9±3)xl0i'' 


Non-local [5] 


OH 


13 


yes 


(1.5-2.5) xlOi'3; 3.2x10^'' 


Non-local [3]; Ground state [4] 


H3O+ 


6 


yes 


1.6x10" 


LTE [6] 


CH 


2 


yes 


9.3x10" 


Ground state [7] 


CH2 


2 


yes 


3.4x10" 


Ground state [7] 


C3 


16 


no 


(4-8)xl0i5 


Non-local [8] 


HF 


1 


no 


1.7x10" 


[9] 


H2D+ 


1 


yes 


9x10" 


[10] 



fl] iGoicoec hea ct ah (2004"); [2] 'Ccc careUi et al. ' 1*2002^: [3] 'Goicocchoa & Cernicharc/ ('2002^; [4] 'Polohampton ct al.' ('2005a'); [5] 
ICernicharo et al. 12006'): [6] Goicocchoa fc Cernicharo, (2001'): [7] Polchampton ct al . (2005b,): [8] Ccrnicharo et a l. (2000): [9] assuming 
the abundance of 3x10"^" from lNeufelT'et'^ 1 1I1997I) and their radial density profile bet ween 0.6 and 22.5 pc. N ote that Neufeld et al. 
used a slightly different beam position on the source to that adopted for the survey; [10] ICernicharo et al.l 1 2007f ). calculated based on 

non-LTE excitation. 



surround ing Sgr B2 using the L WS grating spectrometer. 
However, ICernicharo et al.l (|2006l ) have detected the J=7-6 
Une at 371.65 /xm (806 GHz) from the ground. They per- 
formed non-local radiative transfer models to try to repro- 
duce the lack of CO emission in the LWS range, showing that 
low H2 density is required. If the J=7-6 emission occurs in 
gas with a kinetic temperature ~100 K, the lack of FIR CO 
lines can be reproduced with n(H2) — 2 x 10* cm~^, but if 
the temperature is higher, the density limit is < 10* cm~^. 

The full list of detected lines is shown in Table [S] 
Th is updates an d exte nds t he line lists already pre sented 
bv iPolehamptonI (|2002l ) and ICoicoechea et al.l (|2004l ). The 
remaining unidentified features are described in Sect. 14.71 
The transitions in Table [3] are described by the quantum 
numbers A'^, J and K, where A'^ is the rigid body rotational 
quantum number, J is the total angular momentum quan- 
tum number excluding nuclear spin, and K is the projection 
of the total angular momentum. Further splitting due to nu- 
clear spin is not included in Table O and where this occurs, 
the wavelengths quoted are the average over the split levels 
weighted by the Einstein coefficients. The notation used for 
H2O, CH2 and NH2 is NKaKc, for NH is Nj, for H2D+ is 
JKaKc, and for HsO"*" and NH3 is J^^~ , where +/— speci- 
fies the parity of the inversion state. For C3, the transitions 
are quoted as lines in the P-, Q- and i?-branches of the 1^2 
bending mode. For OH and CH, the lines are quoted for the 
two rotational ladders due to spin-orbit coupling. For the 
low-J levels involved, OH is close to Hund's coupling case 
(a) and the ladders can be labelled ^ni/2 and ^n3/2, where 
the subscript is the total (orbital -I- spin) angular momen- 
tum quantum number of the electrons, Q. CH is very close 
to Hund's case (b) in which there is much stronger coupling 
of electron spin with rotational motion. In this case the elec- 
tron spin along the internuclear axis is not well defined, and 



the ladders are labelled Fi and F2 for the upper and lower 
spin components with a given J-value respectively. The tran- 
sitions for both species are described by the total angular 
momentum quantum number, J. 

Absorption due to the envelope of Sgr B2 occurs at 
a LSR velocity of ~65 km s~^, although other features 
have been observed cen tred at velocities between 50 km s~^ 
and 70 km s"^ (e.g. iMartfn-Pintado et al.l Il990l ). Addi- 
tional absorption occurs in the ground state lines (e.g. for 
OH, CH and CH2) at velocities between —110 km s^^ and 
-1-30 km s~^. This is due to the galactic spiral arms that cross 
the line of sight (e.g. the Galactic Bar, the expanding molec- 
ular ring at 3-4 kpc from the Galactic Centre and the local 
spiral arms; see lCreaves fc Williams 1994 ). Eight individual 
components have been detected in H i absorptio n observa- 
tions towards Sgr B2 M (|Garwood fc Dickevll989l ). although 
at very high spectral resolution, each of these components 
can be separated into many narrow features with velocity 
wi dths ~1 km s~^ (e.g . in th e CS absorption measurements 
of|G reaves fc Williamslll994h . 

In order to determine the characteristics of each de- 
tected feature at the velocity of Sgr B2, we performed a 
fit using the Lorentzian profile of the FP with three free 
parameters describing the line centre, resolving power and 
relative depth below the continuum. The results of these fits 
are shown in Table [3] To improve the accuracy of the fit, 
fine adjustments were made to large scale baseline already 
achieved with the re-binning and spline technique described 
in Sect. 13.41 This was carried out with low order polynomials. 
Note that no correction was made in Table [3] for systematic 
errors in wavelength calibration in the non-prime data be- 
low 70 nm. This can affect the quoted line centres by up to 
25 km s-^ 

This semi-automatic fitting technique was successful 
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Figure 6. Observed rotational transitions of NH3, NH2 and NH. Only the strongest NH3 transitions to metastable levels are shown. 
Note that the central wavelength of the NH3 3^— 2 J transition is shifted with respect to the other lines because it is a blend incorporating 
-2^" transitions. The dashed line shows the velocity of Sgr B2 itself (65 km s~^). 



for the majority of lines. However, some weak features 
either required a more d etailed analysis to fit (e.g. HD; 
IPolehampton et al.l 1200231 ). or L04 data to increase the 
siRnal-to-noise (e.g. OH ' ^'^-1/2 J = 5/2-3/2 at 98.7 /xm; 
iGoicoechea fc Cernicharol |2002| ). These features are indi- 
cated in Table |3] For lines in which line of sight absorption 
was observed, only an estimate for the peak absorption and 
where possible an estimate of the corresponding velocity, is 
reported in the table. More detailed analysis of the line of 
sight features is summarised in the following sections. 

Only a few of the lines have been previously observed 
towards Sgr B2 with lower spectral r esolution using the 
KAO (OH at 119 ^.m; CH at 149 Istorev et all Il98ll : 
IStacev et ai]|l987l ) and so most were observed for the first 
time with ISO. The ISO spectra of the main species have 
already been investigated in detail in previous papers (see 
Tables [5] and 13 for references), however, there are several 
new features whose detection is only reported here. These 
are indicated in bold face in Tabled The following sections 
give an overview of the conclusions for each group of species. 



4.3 Molecular species 

4.3.1 Nitrogen bearing molecules: NHi, NH2 and NH 

The greatest number of lines from a single molecule in the 
survey are due to rotation-inversion transitions of NH3. In 
total there are 21 detected absorption transitions with up- 
per energy levels in the range 65-688 K, including both 
metastable (see Fig. [Sjl and non-metastable lower s tates (see 
detailed results presented bv lCeccarelli et al]l2002l ). Pure in- 
version transitions have been de tected in the radio regio n 
from many of these levels (e.g. iHiittemeister et al.lll993 ). 
recently extending much higher in ener gy t o Jk = 18i8 at 
~3130 K above ground (, Wilson et al.ii2006l ). A simplified 
model of the FIR lines, based on these radio lines showed 
that the higher energy transitions are optic ally thin but the 
lower energy lines have optical depths ~1.6 (|Ceccarelli et al.l 
I2OO2I ). They calculated rotational temperatures of 130±10 K 



for the metastable levels and 310±100 K for the non- 
metastable levels using a standard rotational diagram ap- 
proach. 

Detailed large velocity gradient (LVG) modelling 
showed that the NH3 absorption must originate in a high 
temperature layer in front of the FIR continuum emitted by 
Sgr B2. The best fitting model parameters gave a temper- 
ature of 700±100 K and density < lO'^ cm"^ for this layer. 
Further constraints on the density can be derived from the 
lack o f CO rotational emission in the survey. ICeccarelli et al.l 
l|2002h derive an upper limit on the density of the hot region 
of lO'' cm~^ based on the lack of FIR CO emission, i n goo d 
agreement with that calculated bv lCernicharo et al.l (|2006D . 
The results of the modelling gave a total column density of 
NH3 ~ 3 X 10^^ cm~^, equally shared between the ortho and 
para forms (i.e. a non-LTE ortho/para ratio). 

There is also a clear detection of the closely re- 
lated nitrogen hydrides NH2 and NH in the survey (see 
Fig. [S]). The resu l ts for thes e lines have been present ed by 
ICernicharo et al.l l|2000l ) and IGoicoechea et al.l (|200l). The 
survey shows a detection of both ortho- and para-NH2, the 
first such observation of the ortho form. One para-NH2 line 
(221-I10 J = 1-1) is detected at 126.8 /im. The profile of 
this line is shown in Fig. |6] and is much broader than the 
ortho lines. This is due to a blend with the ground state line 
of ortho-H2D+ (see Sect. BXi)) . 

Para-NH2 has been previously obse rved in Sgr B2 from 
the g round via its mm wave transitions (|van Dishoeck et al.l 
Il993l ). These mm observations indicated a total column den- 
sity of ortho-fpara NH2 of (5.2 ± 1.2) x 10^^ cm"^, assuming 
an ortho/para ratio of 3. It is difficult to use the ISO obser- 
vations to provide measured confirmation of this ortho/para 
ratio due to the uncertainties in the one detected para line. 
However, the fact that the other components of the NH2 
221-I10 transition were not de tected indicates t h at a r atio 
of 3 is a reasonable estimate. IGoicoechea et al] (120041 ) es- 
timated the total column density of NH2, taking account 
of both the ISO and mm data, to be Af(NH2) = (1.5- 
3) X 10^^ cm~^. The column density they estimated for NH 
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Table 3. Identified Lines in the LOS prime data. See text for a description of the notation used to describe each transition. For transitions 
with hne of sight absorption, only an approximate Hne to continuum value is given at the velocity of Sgr B2 (in parenthesis). The line 
to continuum ratio, central velocity and line width were determined by fitting a Lorentzian profile to the L03 survey data (see text). No 
L04 data were used for these fits - features which were only detected in L04 data are noted. References in the literature are shown in 
the last column (both for ISO and KAO detections). References for the line wavelengths are given at the end of the table. Features not 
reported elsewhere arc shown in bold font. 



Rest Species Transition Line v /\v Reference 

Wavelength to (km s"-*^) (km s" -'^) 

(^m) continuum 



51.8145 


[Oiii] 


='P2-^Pi 


1.28 ±0.02 


39 ±3" 


82 ± 7 


[1] 


53.2615 


OH 


''ni/2-2n3/2 J = 3/2+-3/2- 


(0.69) 


~ 47" 


broad 


[2] 


53.3512 


OH 


''ni/2-2n3/2 J = 3/2--3/2+ 


(0.69) 


~ 47" 


broad 


[2] 


56.3387 


p-NHs 


98-8/ 


6 






[3] 


57.3295 


[N in] 


P3/2 — Pl/2 


1.12 ±0.01 


42 ±3" 


50 ±6 




57.6365 


P-H2O 


422 — 3l3 


C, 








62.7274 


P-NH3 




0.960 ± 0.006 


49 ±4" 


41 ± 9 


[3] 


63.1837 


[Oi] 


3Pl-3P2 


(0.50) 




broad 


[4,5,6] 


63.3765 


P-NH3 


87 -7^ 


0.976 ± 0.003 


39 ± 5" 


64 ±15 


[3] 


65.1316 


OH 


203/2 J = 9/2- - 7/2+ 


~ 0.98 


46 ± 6° 


~ 60 




65.2788 


OH 


203/2 J = 9/2+ -7/2- 


0.971 ±0.006 


46 ± 6° 


60 ±20 




66.4377 


0-H2O 


330-221 


0.866 ± 0.01 


49 ± 1" 


45 ± 3 


[7] 


67.0891 


P-H2O 


331-220 


C 






[7] 


69.5377 


H3O+ 


4J-33 


0.950 ±0.03 


51 ± 14" 


36 ±28 




71 finSA 

1 J-.UU04: 


0-NH3 


7"^ fi~ 
'6 ^6 


u.ouu in u.uz 


4- 9 
00 in z 


97 ^ ^ 






O-iN XI3 


7~ F,~^ 

'6 


n Qni 4- n 09 


f^S 4- 9 
00 III Z 






/ Z.oZoo 


p-lN J:l3 


"5 -65 


b 






[3] 


( O.ooU ( 


0-rl2 W 


021— '212 


u.Dyo It u.uuy 


c^O 4- 1 
O J It 1 


/I Q 4- 9 
4o It Z 


^7^ 
I'i 


f o. / 


, . II, r\ 


423 — 021 


n 4- n ni 


p;q 4- 
oy It 


AQ 4-11 

44y It X X 




VO 1 1 IPi 

/ y. ii (0 






(O.b) 


~ 00 


broad 




7Q 1 9 


ow 

wn 


2n /„ 2tt„,„ 7 — 1/9- 






broad 


[91 


82.2742 


p — H3O + 


^2 ~ -^2 


(0.875) 


~ 57 


broad 




83.4320 


P-NH3 




0.860 ± 0.01 


59 ± 6 


42 ± 9 


[3] 


83.5898 


P-NH3 


6J-54 








[3] 


84.4201 


OH 


2n3/2 J = 7/2+-5/2- 


0.801 ± 0.01 


58 ± 1 


37 ± 4 


[2,8] 


84.5441 


P-NH3 




0.885 ± 0.04 


58 ± 5 


26 ± 8 


[3] 


84.5966 


OH 


2n3/2 J = 7/2--5/2+ 


0.802 ± 0.03 


53 ±4 


37 ± 7 


[2] 


88.3564 


[Om] 


^'Pi-^'Po 


1.193 ± 0.01 


69 ±1 


52 ± 2 


[1] 


89.9884 


P-H2O 


322-2ll 


0.897 ± 0.03 


52 ± 3 


36 ± 5 


[7] 


98.7310 


OH 


^ni/2 J = 5/2-3/2 


d 






[2] 


99.4931 


- H2O 


5o5 — 4i4 


0.987 ± 0.005 


46 ± 4 


46 ± 10 




99.9498 


P-NH3 


5t-47 
4 4 


0.811 ± 0.007 


57 ± 1 


37 ± 2 


[3] 


100.1048 


0-NH3 




0.962 ± 0.03 


54 ±2 


49 ± 9 


[3] 


100.2128 


P-NH3 


5+-4- 
■^2 ^2 


b 






[3] 


100.5767 


P-H3O+ 


2+-1- 


(0.905) 


~ 60 


broad 


[9] 


100.8686 


0-H3O+ 


2+-1- 
^0 -^0 


(0.860) 


~ 60 


broad 


[9] 


100.9828 


P-H2O 


220-111 


0.695 ± 0.006 


64 ± 1 


51 ± 1 


[7,9] 


101.5337 


P-NH3 


54" -4 


0.819 ±0.005 


60 ±1 


42 ± 2 


[3] 


101.5965 


0-NH3 


53" -4 


0.976 ± 0.004 


57 ±3 


31 ± 7 


[3] 


102.0050 


HfO 


220-111 


0.963 ± 0.01 


49 ±3 


32 ± 9 


[7] 


107.7203 


P-CH2 


2i2-ll0 J = 3-2 


(0.97) 




broad 


[10] 


108.0732 


0-H2O 


221-I10 


0.622 ± 0.04 


63 ±1 


49 ± 1 


[7] 


109.3466 


HiSQ 


221 — lio 


0.956 ±0.008 


61 ± 5 


40 ±12 




112.0725 


HD 


J = 1-0 


e 






[11] 


113.5374 


- H2O 


4l4 — 3o3 


0.944 ±0.007 


58 ±2 


31 ±6 




117.0648 


0-NH2 


220-I11 J = 2-1 


0.936 ±0.006 


63 ±2 


45 ± 7 


[1] 


117.3831 


0-NH2 


220-I11 J = 2-2 


0.889 ± 0.006 


63 ±1 


49 ±4 


[1] 


117.7918 


0-NH2 


220-I11 J = 1-1 


0.864 ±0.01 


64 ±2 


38 ±4 


[1] 


119.2334 


OH 


2n3/2 J = 5/2-3/2 


(0.14) 




broad 


[2,8,12] 


119.4409 


OH 


2n3/2 J = 5/2-3/2 


(0.14) 




broad 


[2,8,12] 


119.621 


i^OH 


2n3/2 J = 5/2-3/2 


(0.97) 




broad 


[13] 


119.828 


i^OH 


2n3/2 J = 5/2-3/2 


(0.97) 




broad 


[13] 



" note that no correction for systematic errors in wavelength calibration for non-prime data below 70 fim has been applied 
''too weak to fit successfully, but observed by ^eccarell^^t_^L_ (2002) ; fit unsuccessful 
too weak to fit successfully, but observed in L04 data by ^oicoeche^_et_aL (2004); see £oleham£toi^_et_^l^ (2002a) 
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Table 4. Identified Lines in the LOS prime data continued. 



Rest 


Species 


Transition 


Line 


V 




Reference 


Wavelength 






to 


(km s ) 


(km s~l) 




{lira) 






continuum 








119.9664 


18 OH 


2n3/2 J = 5/2-3/2 


(0.9) 


~ 65 


broad 


[2,13,14] 


120.1730 


18 OH 


^n^/o J = 5/2-3/2 


(0.9) 


~ 61 


broad 


[2,13,14] 


121.6973 


HF 


J = 2-1 


0.971 ± 0.006 


66 ± 5 


28 ± 8 


[15] 


121.8976 


[N ii] 


_3p^ 


/ 








124.6475 


0-NH3 


4+-3,7 

3 A 


0.755 ± 0.004 


61 ± 1 


47 ± 1 


[31 


124.7958 




^2 "^2 


0.973 ± 0.004 


57 ± 4 


33 ± 8 


f3l 
L"J 


124.8834 


n-NHc 


4+-3- 
^1 '-'1 


951 + 01 


57 ± 3 


21 ± 11 


[31 


1 26 801 4 


P-NH2 


2oi — 1 1 n 7—1-1 


qf!c; + n 01 


64 ± 3 


68 ± 19 


111 


126.8530 


0-H2D+ 


2l2 111 


blend 




blend 


fl9l 


127.1081 


o-NHq 


4--3+ 
^3 '^3 


Q 7QQ + f) 01 


68 ± 1 


34 ± 4 


[31 






1 1 1 C[nn 7 — 9 1 

ril UOO 'J — ^ ^ 






bro&d 


[ml 


1 27 8582 


0-CH2 


111 -Onn 7—1-1 


(0.98) 




brocid 


fiol 




^3 






64 ± 6 


42 di 14 




131.9629 




R (12) 


0.961 di 0.008 


65 ± 4 


56 di 14 




XOt-J ■ C 7 1 J \_J 




R(10) 




67 ± 4 


31 di 12 




1 38 5278 


n-H-iO 
P 1^2 ^ 


3i3~2q2 


n 031 + n n08 


64 ± 2 


25 di 5 


171 
L'J 








9 








143 8801 




R(6) 


941 + 003 


69 ± 1 


43 di 4 


1161 






^0 ^1 


1 1 74 _|_ n nofi 


75 ± 1 


54 di 3 


11 fil 


1 48 041 9 


^3 




n Q1 f; _|_ n nny 

L/.jJ-TJ _l_ yj.yjyji 


70 ± 2 


39 di 5 


[lOJ 


149.0912 


CH 


Fo 3/2-1 /2 


(0 641 


~ 65 


brocid 


[1,10,17] 


149 38Q9 


CH 




fO 64") 


~ 65 




[1,10,17] 


151 5279 


NH 


2i — In 


0.962 ± 0.01 


58 ± 10 


43 ± 24 




152.2875 


C3 


R(2) 


0.943 ± 0.01 


84 ± 3 


25 ± 7 


[16,18] 


153.0961 


NH 


22-I1 


0.925 ±0.02 


69 ± 4 


31 ± 11 


[1] 

L J 


153.2956 


C3 


Q{12) 


s 






[16] 


153.3444 


NH 


23-I2 


s 






[1,16] 


154.8619 


C3 


Q{10) 


0.973 ±0.01 


55 ±7 


45± 15 


[16] 


156.1898 


C3 - blend 


Q(8) 


0.939 ±0.01 


69 ±3 


28 ±7 


[16] 


156.1940 


P-H2O - blend 


322-3l3 


blend 


(61 ± 3) 


blend 


[16] 


157.2609 


C3 


Q(6) 


0.878 ±0.01 


67 ±3 


52 ±67 


[16] 


157.7409 


[Cii] 


2p , 2p 
-^3/2^ ^1/2 


1.655 ±0.01 


74 ± 1 


49 ± 1 


[1,6] 


158.0595 


C3 


Q(4) 


0.894 ±0.01 


72 ± 3 


33 ±6 


[16] 


158.5735 


C3 


Q(2) 


0.937±0.02 


68 ±3 


~ 18 




163.1247 


OH 


2Hi/2 3/2-1/2 


1.065 ±0.01 


65 ±4 


41 ± 11 


[2] 


163.3973 


OH 


2Hi/2 3/2-1/2 


1.060 ±0.01 


75 ±4 


35± 10 


[2] 


165.5966 


p-NHa 


■^2 ^^2 


0.711 ± 0.01 


64 ± 1 


54 ±4 


[3] 


165.7287 


P-NH3 


3+-2- 


0.911 ±0.03 


82 ±4 


44± 13 


[3] 


167.6794 


C3 


P(4) 


0.948 ±0.01 


70 ±7 


54± 15 


[16] 


169.9676 


P-NH3 - blend 


•^2 ^2 


blend 


85 ± 1 


blend 


[3] 


169.9887 


P-NH3 - blend 


3r-2+ 


0.755 ±0.01 


48 ± 1 


62 ±4 


[3] 


169.9961 


0-NH3 - blend 




blend 


35 ± 1 


blend 


[3] 


174.6259 


0-H2O 


3o3-2l2 


0.770 ±0.01 


65 ± 1 


32 ±2 


[7,8] 


179.5267 


0-H2O 


2l2-loi 


(0.08) 




broad 


[7,8] 


180.2088 


0-H3O+ 


■^3 - ^3 


(0.84) 


~ 60 


broad 




180.4883 


0-H2O 


221-2i2 


0.715 ±0.03 


63 ± 3 


48 ±9 


[7] 


181.0487 


0-HI80 - blend 


2l2-loi 


0.714 ±0.03 


64 ±2 


36 ±7 


[7,8,9] 


181.0545 


P-H3O+ - blend 




blend 


(55 ± 2) 


blend 


[7,8,9] 



i broadened due to hyperfinc structure; ^ contained within broader absorption profile 

Line wavelengths: [Olll]: Pcttc rsson f 1982 1: OH: Brown et al. ('1982' ): Varbcrg & Evcnson (1993); NH3: JPL catalogue; 

[Nm]: Gry ct al. (2003); [Ol]: ZinketaL (1991); H9O: Matsushima et al. (1995): Johns (1985); Hi^O: Matsushima et al. (1999) : 

H3O+: JPL catalogue; CH9: Polehampton et al. (2005b'): HP: Evenson et al. (1988) : NH7: Morino fc Kawaguchi (1997) : 

I'^OH: Polehamp ton et al. (2003) : 18 OH: Morino et al. (1995) : HF: Nolt et al. (1987) : [Nil]: Brown et al. (1994) ; C3: Giesen et al. (2001) ; 

CH: Davidson ct aL (2001) : NH: JPL catalogue; fCll]: Cooksv et al. (1986): H9D+: Cernicharo et al. f2007) . 

References: [1] Goicoec hea et al. (2004) ; [2] Goicoechea fc Cernicharo (2002) : [3] Ceccarelli et al. (2002) ; [4] Baluteau et al. (1997) ; 

[51 Lis et al. (2001) : [6] Vas tel et al. (2002) : [71 Cernicharo et al. (2006) : [81 Cernicharo et al . (1997a); [9] Goicoechea fc Cernicharo f2001) : 

[10] Polehampton et al. (2005b) : [11] Polehampton et al. (2002a) ; [12] Storev et al. (1981) ; [13] Polehampton et al. (2003) : 

[14] Lueten et al. (1986) ; [15] Neufeld et al. (1997) : [16] Cernicharo et al. (2000) ; [17] Stacev et al. (1987) ; [18] Giesen et al. (2001) ; 

[19] Cernicharo et al. (2007) . 
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Figure 7. Observed rotational transitions of para and ortho water. The energy level diagram for water is shown in the centre, with 
transition wavelengths given in microns. Not all of the para lines shown were detected. The dashed line shows the velocity of Sgr B2 itself 
(65 km s~^). The insert on the right is the ground state ortho-H20 line, which shows complete absorption at all line of sight velocities. 



was ~ 4 X 10^'* cm~^, giving a final ratio NH3/NH2/NH ap- 
proximately equal to 100/10/1. Goicoechea et al. show that 
these r atios cannot be ex plained by either dark cloud mod- 
els Ce.g.lMillar et al.lllQQll'). o r by UV illuminated PDRs (e.g. 
ISternberg fc Dalgarno|]l995h . However, low- velocity shocks 
can heat the gas up to 500 K (consistent with the tempera- 
tures derived from ISO NH3 observations) as well as satisfac- 
torily reproducing the observed NH3 and NH2 abundances, 
providing that ammonia is efficiently formed on grain sur- 
faces and sputtered back t o the gas phase by t he shock pas- 
sage (see detailed models of lFlower et al.|[l995l ). This implies 
that shocks contribute to the energy balance of the Sgr B2 
envelope as well as other mechanisms such as those associ- 
ated with the presence of far-UV radiation fields. 

All three nitrogen species only show absorption at 
the velocity of Sgr B2, and no line of sight absorption 
is observed. However, the ground state ammonia line has 
been observed by the ODIN satellite towards Sgr B2 
marson et alj|2005l ). This does show absorption along 
the whole line of sight. 



4--3.2 Oxygen bearing molecules: H2O, OH and HsO'^ 

Water is an extremely important species in the survey as it 
has the largest associated column density of all the detected 
molecules. It is also an important research field in its own 
right, particularly in the warm neutral gas driven by star 
formation, where it plays a dominant role in the thermal 
balance. As water is an abundant species in the Earth's at- 
mosphere, only space telescopes have access to its thermal 
rotational lines and only very limited studies were possible 
before ISO. 

Many rotational lines of water have been observed in the 
survey. These are mostly due to Sgr B2 itself, although in the 
lowest energy lines there is some absorption at the velocities 
of line of sight clouds. The water transitions detected in the 
survey are shown in Fig. [7] These lines ha ve been analysed 
in detail bv lCernicharo et al. I l|l997al . l2006D using data from 
the LWS L04 mode. They presented data for 14 water lines 



and 2 lines of H2*0 observed with ISO, and combined these 
with a map of emission by the 813-220 183.31 GHz maser 
line observed with the IRAM 30-m telescope. The FIR lines 
are all in absorption and optically thick (particularly the 
2i2-loi at 179.5 /xm which has an optical depth of ~ 10'^- 
10*). The H2O absorption traces the outer surface layers of 
the Sgr B2 envelope (i.e. the same hot, low density layer 
observed in NH3), whereas the 183 GHz line traces denser 
gas closer to the hot cores in the inside of the cloud. 

ICernicharo et af] l|2006l ) carried out LVG and non-local 
radiative transfer modelling and found that IR photons from 
the dust are the dominant source of excitation of the water 
rotational levels. This showed that the observed lines were 
not very sensitive to the gas temperature, and detection of 
weaker lines that would be more sensitive are needed to bet- 
ter constrain the physical conditions of the water absorbing 
layers. However, assuming that t he FIR water lines arise in 
the s ame warm gas traced by OH (|Goicoechea Sz Cernicharol 
|2002| ). the best fitting column density for the ISO observa- 
tions was (9 ± 3) X 10^** cm^'^, with an abundance of (1- 
2)xl0~^ 

Previous observations of the ground state water line 
with SWAS, and HDO lines with ground ba sed telescopes 
(iNeufeld et al. I l2000l . I2OO3I : IComito et ai]|2003D indicate that 
there is also a component of the absorption that is due to 
the warm envelope surrounding the Sgr B2 cores as well as 
the hot layer observed in NH3 lines. If the majority of the 
water absorption is located in the hot foreground layer, the 
estimated column density of H2O associated with Sgr B2 cal- 
culated from the SW AS observations is (2.5-4) x 10^® cm"^ 
(|Neufeld et al .1120031 ) . of the same order of magnitude as that 
calculated from ISO bv lCernicharo et al. 

Cernicharo et al. I (Il997al). iGoicoechea Sz Cernicharol 



(|2002t ) and I Goicoechea fc Cernichard (|200li r calculate an 
abundance for H2O of a few 10~®. This abundance gives 
a ratio with OH of H20/OH=2-4. This relatively low ra- 
tio can be reproduced if the Sgr B2 envelope is illumi- 
nated by a strong far-UV radiation field, and the gas is 
warm enough to activate neutral-neutral reactions involv- 
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Figure 8. Observed rotational transitions of OH within and be- 
tween tile s-iid ^113/2 ladders. The ground state line in the 
■^113/2 ladder is shown for both ^^OH and ^**OH. The dashed line 
shows the velocity of Sgr B2 itself (65 km s~^). 



ing H2O and OH. The presence of such a radiation field is 
inferred from the extended atomic and ionic fine structure 
line emission (see Sect. 14.4] ). Detailed photochemical models 
adapted to the conditions of the Sgr B2 envelope reproduce 
the above ratios if most of the observed FIR H2O and OH 
absorption arises from t he photoactive surface of the cloud 
ICernicharo et ahl I2OO6I ) . Note that the same shock mod- 
els that apparently reproduce the observed NH3 and NH2 
abundance ratios (Sect. 14.3. ip fail to reproduce the results 
for water and related species (at least in the averaged picture 
contained within the large ISO beam). In particular, the pre- 
dicted H2O column density is almost 2 orders of magnitude 
larger than observed by ISO, while the predicted OH col- 
umn density is an order of magnitude lower than observed. 
In the line of sight clouds, the ratio of H 2O to OH is even 
lower, 0.6-1.2 jPolehampton et al.ll2005al ). almost as low as 
the ratio calculated for line of si ght features towards W51 
and W49 of 0.3-0.4 l|NeufeId et a l. 2002 ; Plume ct al. 2004). 

The survey also contains lines due to the pure rotational 
transitions of the OH molecule up to 400 K above ground 
(see Fig. [H]). Only the strongest ground state transition had 
previously been ob served towards Sgr B2, using the KAO 
IStorev et al.|[l98ll ). The ground state lines show absorption 
due to the entire line of sight, but at higher energies only 
Sgr B2 itself is detected. The detected OH transitions have 2 
components due to the A-doublet splitting of each rotational 
level, and in most cases these were resolved by the LWS FP. 
No asymmetry in the absorption was observed between any 
of t hese doublets. 

iGoicoechea fc Cernicharcl (|2002l ) have carried out a de- 
tailed radiative transfer analysis of these lines, and conclude 
that they also arise from the low density (<10^ cm~'^) en- 
velope of Sgr B2 where the temperature should range from 
Tfc=40 K in the innermost regions to Tfc=600 K in the cloud 
edge. The high abundance of OH, ~(2-5)xl0~^, and large 



OH/H2O abundance ratio is attributed to the presence of 
clumpy PDRs on the edge of the Sgr B2 envelope. 

Here, we present a weak detection of an additional ab- 
sorption doublet of OH a t 65 ^m (^n3/2 J=9/2-7/ 2 ). Ac - 
cording to the models of IGoicoechea fc Cernicharol ( 20021 ) 
this line is only predicted to be in absorption if it arises from 
the colder regions of the envelope, or if the FIR continuum 
level at ~65 /im was underestimated for the warm gas. The 
wide range of conditions traced by the observed OH lines 
mean that the additional 65 /^m detection does not provide 
significant new constraints on the previous models. 

Overall, the OH excitation is dominated by FIR pump- 
ing, with transitions in the ladder seen in emission. 
This is due to the upper levels being populated by the cross 
ladder transitions at 53 /im (^ni/2-'^n3/2 J=3/2-3/2) and 
79 ^im (^ni/2-^n3/2 J=l/2-3/2) - see Fig.U 

In addition to absorption from Sgr B2, the transitions 
from the ground rotational state show absorption at veloc- 
ities associated with clouds along the line of sight. This 
includes the two cross ladder transitions at 53 fim and 
79 fim and the fundamental transition in the ^n3/2 ladder at 
119 nm, which is highly optically thick at all line of sight ve- 
locities. The higher energy transitions are not seen in these 
clouds, indicating that they are colder with lower excita- 
tion. The line of sight components were fitted in the 53 fj.m 
and 79 /^m lines by constructing a high resolution model of 
10 velocity components, based on the line widths and ve- 
l ocitie s from the H I measurements of [Garwood fc DickevI 
(Il989l ). This model was then convolved to the resolution of 
the LWS FP (see IPolehampton et all bOOSal ). As only the 
ground state lines were observed in the line of sight clouds, 
the total column density of OH can be estimated simply from 
the best fitting high resolution model. The optical depth for 
each component, r, was adjusted to find the best fit with 
the line to continuum ratio set by. 



I = /cBxp (-r) 



(1) 



where Ic is the intensity of the continuum. Ground state col- 
umn densities for each component (in cm~^) were calculated 
assuming a Dopplor line profile with Maxwellian velocity dis- 
tribution (Spitzcr 1978), 



tAv 



(2) 



Aij is the Einstein 



where Av is the line width in km s~^ 

coefficient for spontaneous emission in s~^, Xij is the wave- 
length in /xm and gi is the statistical weight of state i. The 
fitted components were all optically thin, except for the one 
at the velocity of Sgr B2. The column density for each line 
of sight feature are shown in Table |S] The OH abundance 
in the line of sight clouds was estimated using molecular 
hydrogen column densities from .Gre aves fc Nyman ( 1996|) 
to be in the range 10~^- 10~^, whereas at the velocity o f 
Sgr B2, it is (2-5) x 10"" (jOoicoechea fc Gernicharoll2002l ). 
OH absorption also occurs across the extended region sur- 
rou nding Sgr B2 as seen b y the LWS grating spectrome- 
ter (jGoicoechea et al] 120041'). Th e final column density de- 
termined by Polehampton et"al] l|2005ah at the velocity of 



Sgr B2 was 3.2 



+0 '^ X lO^"^ cm-2 



in reasonable agreement 



with iGoicoechea fc Cernicharol (|2002l ) who derived a value 
of (1.5-2.5) X 10^® cm-^ 



The ISO LWS spectral survey towards Sgr B2 13 




Figure 9. Observed rotation-inversion transitions of HsO^. The energy level diagram in the centre shows transition wavelengths in 
microns. Note that the para line is blended with Hj^O. The dashed line shows the velocity of Sgr B2 itself (65 km s~^). 



Both the isotopic species of OH, ^'^OH and ^*0H, were 
also detected in the survey and are described in Sect. 14.61 

Another species closely linked to water and hydroxyl 
is H3O'''. This molecule is similar in structure to ammonia 
but has a very large inversion barrier, meaning that it shows 
pure inversion transitions in the FIR. The results for several 
of the lower energy lines have already been presented by 
iGoicoechea fc Cernicharol (|200ll ). and show absorption due 
to the entire line of sight. Here, we present several additional 
detections of H3O''' lines at shorter wavelengths and higher 
energy levels. The detection of higher excitation HsO"^ lines 
in absorption indicates that HsO^ is also likely to arise in 
the warm gas traced by H2O and OH in the Sgr B2 envelope. 
Further detailed analysis and modelling is required. 

Both water and OH are a product o f HaO"*" dissocia- 
tive re combination. However, as noted by ICernicharo et al.l 
(|2006l '). the large H2O and OH column densities found to- 
wards Sgr B2 can only be reproduced if a significant frac- 
tion of the gas is warm enough {Tu ~300-500 K) to acti- 
vate additional neutral-neutral reactions that can efficiently 
form larger amounts of water and OH than HsO^ disso- 
ciative recombination alone. Detailed photochemical models 
assuming that the Sgr B2 envelope is illum inated by 10^-10^ 
times the mean interstellar ra diation field (jGoicoechea et ahl 
I2OO4I : ICernicharo et all |2006| ) satisfactorily reproduce the 
H2O/OH abundance ratios and absolute column densities, 
as well as the Ol line intensities (see Sect. |4^ . Therefore, 
the oxygen chemistry seems to be dominated by the pres- 
ence of warm PDRs in the external layers of Sgr B2 (at 
least in the averaged picture provided by the large beam of 
ISO observations). 

The region around Sgr B2 has recently been mapped in 
the 3^-22 (364 GHz) and 1 7-2+ (307 GHz) lines o f H3O+ 
using the APEX telescope (|van der Tak et al.ll2006l '). How- 
ever, these lines appear in emission and trace denser gas than 
the absorption lines presented here in the core of Sgr B2. The 
line ratio indicates a high excitation temperature, and at the 
conditions of the core of Sgr B2 M give a total H3O+ col- 



umn density of 1.4 x 10^® cm ^, much higher than the value 
estimated for t he envelope of Sgr B2 from the ISO lines: 
1.6 xlO^"* cm"^ (|Goicoechea fc Cernicharol I2OO h i. 



^.3.3 Carbon bearing molecules: CH, CH2 and C3 

There is a clear detection in the survey of the two A-doublet 
components of the J = 3/2-1/2 transition of CH from 
the ground to first rotational state of the F2 ladder, at 
149 ^m. These line s were first observed using the KAO by 
IStacev et af] l| 19871 ). The results from t he ISO survey have 
been presented by iPolehampton et al.l (l2005bl) and in L04 
observations by ICoicoechea et al.l ( 20M ). Goicoechea et al. 
also observed CH in the extended Sgr B2 region using the 
LWS grating mode. 

In contrast to OH, no higher energy rotational tran- 
sitions of CH have been detected. The next highest energy 
transition in the F2 ladder, ~96 K above ground, is J = 5/2- 
3/2 at 115 /im. Taking an uppe r limit for the absorption i n 
this hne shows that r^t < 20 K jPolehampton et al.|[2"005bh . 
The lowest energy line in the F\ ladder of CH at 180 fim 
(J=5/2-3/2) is also not observed. The lower energy level of 
this transition would be populated by absorption in the cross 
ladder trans ition at 560 /xm (-F1-J2 J=3/2 -l/2). However, 
as noted by IConzalez- Alfonso et al] ((2OOJ), the continuum 
in Sgr B2 at 560 /xm is not strong enough to produce a 
significant population in the F\ ladder (in contrast to the 
ultraluminous galaxy Arp220). 

The lack of absorption from higher energy levels indi- 
cates that the ground state transition should give a good 
measure of the total column density. A similar method 
to that applied to OH was used by IPolehampton et al.l 
(2005b) to determine the column density of CH in the line 
of sight features in the range (0.3-3.1) x 10^'* cm~^, and 
(9.3±0.9)xl0" cm"^ i n the Sgr B2 component . The results 
are shown in Table (6] IGoicoechea et al. 1 l|2004h found inte- 
grated column densities of (0.8-1.8) x 10^^ cm~^ across the 
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Figure 10. The observed P, Q and R branch Unes of C3 in the survey. Data are shown around a wide range of hnes, not all of which 
are claimed as detections. The nominal velocity of Sgr B2 (65 km s~^) is shown by the dashed line. 



extended region, peaking at the central Sgr B2 M and Sgr B2 
N positions. 



The survey also incudes the lowest energy rotational 
transitions of CH2, the det ection of which was described by 
iPolehampton et al.l l|2005bf ). This represents the first detec- 
tion of these low energy transitions and is only t he sec ond 
definitive detection of CH2 in space: iHolhs et all l|l995h de- 
tected emission from CH[2 towards Orion KL and W51M 
via its 4o4-3i3 transition at 68-71 GIfz. Other than this, 
only a tentative assignment of CH[2 absorption bands in the 
UV spe ctrum towards H P 154368 and C, Oph has been re- 
ported JLvu et al.ll200ll ). Both CH and CH2 show broad 
absorption due to the entire line of sight towards Sgr B2. 
They are chemically related as both are formed by the dis- 
sociative recombination of CH|j" and destroyed by reaction 
with atomic oxygen. The shape of the CH2 line could be 
well fitted by assuming a constant CH/CH2 ratio along the 
whole line of sight and the best fit gave CH/CH2 equal to 
2.7±0.5, contrary to the expectation of the CH^ branching 
ratios measured in the lab (C H: 30%, CHj: 40%, C: 30%; 
IVeibv-Christensen et al.lll997l ). In order to explain the low 
abundance of CH2, other formation/destruction routes are 
required. One solution is to include a high UV radiation 
field (possibly applicable in Sgr B2 itself), or to include the 
formation of hydrocarbons on dust grain surfaces (e.g. the 



models of lViti et al]|200ll . give CH/CH2 ratios of 7.7-13.6 
for Av=4). 

Finally, 16 lines of triatomic carbon, C3, have been 
detected in the survey. C3 has its lowest energy bending 
mode (i^2=l-0) centred w ell within the ISO su rvey range 
at 157.69 /xm (63.42 cm~^: [Schmuttenmaer et al.i,1990 ) and 
this is the only molecule in the survey detected through vi- 
brational transitions. 

C3 has previously been observed in the C-rich evolved 
star IR C-l-10216 through i ts high energy stretching mode at 
4.9 (im. (|Hmkle et al .1 19881 ). as well as via emission in some of 
the s ame FIR transitions as seen in Sgr B2 (I Cernicharo et ahl 
|2000D . C3 has also been observed in optical absorption spec- 
tra of diffuse and transluc ent clouds in the sight lines of stars 
(|Adamkovics et al]l2003l . and references therein). However, 
in molecular clouds where the expected mid-IR flux is too 
low to allow systematic searches for the 4.9 pim stretching 
mode, and the optical extinction is very high, the FIR bend- 
ing mode is the best way to detect C3 and other such nonpo- 
lar species (e.g. the Cn carbon chains). These species do not 
have a permanent di pole moment and so no pure rotational 
lines. As noted by Ce rnicharo et al] (|2000(). the low- lying 
vibrational bending modes of polyatomic molecules could 
dominate the spectra to be observed by future FIR space 
missions, which will have higher spectral and spatial resolu- 
tion than ISO. 
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C3 has been analysed in detail by ICernicharo et all 

l|2000l '). who observed 9 of the FIR lines using the LWS L04 
mode. They used a non-local radiative transfer model to 
calculate a fractional abundance of ^(Ca) = 3 x 10"* (or 
even higher if C3 only arises from the warm surface of the 
cloud) . In the LOS survey, we detect several more higher en- 
ergy transitions of C3, bringing the total number of lines 
to 16. Figure [10] shows the survey data centred on the C3 
wavelengt hs from the Cologn e Database for Molecular Spec- 
troscopy l|Miiller et al.ll2005l ) up to P(8), Q(16) and i?(18). 
Only 3 lines in the P-branch are detected, but in the Q- 
and i?-branches some absorption is observed in lines up to 
Q(14) and _R(16) (Fig.[TO]also shows the data around several 
additional transitions which are not detected). 

The R{2) line of C3 has also been observed usin g the 
KAO with high spectral resolution l|Giesen et al.|[200ll ). This 
line appears narrow (FWHM of 8.3 km s~M and centred 
at the velocity of Sgr B2 at 63.7 km s"^. ICernicharo et af] 
l|200(t ) indicate that the R{2) line at 152.3 /im may show 
absorption due to the entire line of sight, whereas the other 
lines are only detected at the velocity of Sgr B2 itself. No 
broad absorption is detected in the L03 survey data pre- 
sented here for any of the detected C3 features. The discrep- 
ancy can be explained by the small wavelength coverage in 
the L04 data meaning that it is difficult to determine the 
true continuum level. The narrow absorption is also con- 
firmed by the KAO observation of th is line, which show s 
only the Sgr B2 velocity component (|Giesen et al.l [2OO1I ). 
However, the central velocity observed by the KAO was 
63.7 km s~^, which is in approximate agreement with the 
other C3 lines observed in our survey, but does not match 
our R{2) line velocity. This may indicate that the assign- 
ment for the R{2) line in the survey may be wrong. The line 
observed by ISO occurs centred at a velocity of 84 km s"^, 
higher than all the other observed lines. This could either be 
due to blending with a line from another (unknown) species, 
or due to some contribution from a spurious instrumental ef- 
fect very close to the expected C3 wavelength. 

^.3.4 Hydrogen Fluoride 

The survey shows a clear detection of the J=2-l rotational 
transition of HF at 121.6973 ^m. This line was also ob- 
served with higher s ignal-to-noise using the LWS L04 mode 
l|Neufeld et al.lll997^ . with the telescope pointed half way 
between Sgr B2 M and Sgr B2 N (the L03 survey data were 
pointed so as to only include Sgr B2 M in the beam). The 
absorption line in the survey agrees very well in velocity with 
the L04 data (66 ±5 km s~^ compared with 67 km s~^) but 
the equivalent width is lower (0.5 ± 0.2 nm compared with 
1.0 nm). This difference in equivalent width could indicate 
that the absorption varies with the position of the beam, re- 
lated to difference between the M and N sources. There may 
also be a contribution to the line from the ortho-H20 432-423 
trans ition at 121.71 9 1 ^J,r^^ (separation of ~ 54 km s~^ from 
HF). iNeufeld et all (Il997l) attribute a 5a emission feature 
to this line. However. ICernicharo et al] (|2006l ) have fitted it 
as an absorption component. In any case, the contribution 
must be weak because the peak absorption occurs centred 
at the expected HF position rather than shifted towards the 
H2O wavelength and the signal-to-noise in the L03 survey is 
not enough to distinguish any H2O contribution. 



Table 5. Summary of atomic species observed in the survey. The 
survey fluxes were calculated using the line to continuum ratio 
from Table [3] with an estimate of the continuum flux from the 
LWS grating mode observation. The last colu mn shows the line 
fluxes measured using the LWS grating mode bv lGoicoechea et al.l 
120041'). 



Species 


Estimated 


Survey Flux 


Grating Flux" 




Continuum 


10-18 


10-18 




10-15 


W cm-2 


W cm-2 




W cm"^ 






Om 51.815 


0.66±0.07 


4.1±0.6 


2.2±0.6 


Nm 57.330 


0.81±0.08 


1.5±0.3 


<1.7 


Oi 63.184 


l.OitO.l 


-absn- 


-absn- 


Om 88.356 


1.3±0.1 


6.0±0.6 


4.2±1.7 


Nn 121.898 


0.77±0.15 


3.7±0.7 


<4.4 


Oi 145.525 


0.54±0.11 


3.9±0.8 


7.3±0.8 


Cii 157.741 


0.49±0.1 


13.0±2.7 


b 



" from lGoicoechea et al] 12004') 

' not detected in the grating due to combination of absorption 
and emission 



iNeufeld et all (|l997l ) used an excitation model to deter- 
mine the HF level populations as a function of position in the 
source. HF has a high critical density for excitation of the 
J=l level by collisions and so radiative pumping dominates. 
They assumed a radial density and temperature profile for 
the envelope of Sgr B2 to calculate an HF abundance of 
3 X 10"^". They show that HF is expected to be the major 
reservoir of fiuorine and so this result indicates a depletion 
of fluorin e by a factor of ^50 fr om the Solar System value. 
However. ICeccarelli et all l|2002l ) suggest that there could be 
a contribution to the absorption from the hot layer seen in 
NH3. This would lead to a higher abundance in that layer 
and so a lower depletion factor (althou gh this would increa se 
the depletion in the core of the cloud; INeufeld et al]|2003D . 

The related molecule HCl, has been observed using 
the KAO via its J = 1—0 transition at 478.96 /xm 
dZmuidzinas et cil]|l995l ). However, the lowest energy tran- 
sition of HCl in o ur survey range is J = 3 — 2 at 159.785 
(|Noh et al.lll987l ) and no feature was detected at this wave- 
length. 



4.4 Atomic and ionic lines 

The survey range includes the important atomic cooling 
lines from O I and C 11, as well as lines from ionised oxy- 
gen and nitrogen. Figure [TT] shows the observed atomic 
lines on an absolute flux scale. The continuum flux was 
determined from the LWS grating observation of Sgr B2. 
The atomic lines were also observed with the LWS grat- 
ing in the extended region surrounding S gr B2 where they 
show widespread emission over 9' x 27' jCoicoechea et all 
l2004h . However, the higher resolution FP observations are 
very important for the central Sgr B2 M position where the 
weak lines such as Nlll (57.317 /im) and Nil (121.898 /tm) 
were not detected at the spectral resolution of the grating. 
Also, the lines of Ol (63.184 /^m) and C 11 (157.741 /im) 
show structure due to absorption along the line of sight 
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IVastel et al.ll2002l ) which was not resolved with the LWS 
grating. The fitted emission line fluxes are shown in Table [S] 



4.4.1 Ol and Cll 

The survey spectral range contains the 3 important cool- 
ing lines that trace photodissociation regions (PDRs) at 
the interfaces of molecular clouds, O I ^Pi-^P2 63.2 /^m, 
Oi ^Pq-^Pi 145 . 3 nm and C 11 ^P3/2-^Pi/2 157.7 /^m. 
iGoicoechea et ahl (|2004l ) noted that the shocked gas in the 
Sgr B2 envelope only makes a minor contribution to the 
observed lines fluxes, and therefore the extended O I and 
C II emission is dominated by the PDR scenario. In partic- 
ular, comparison of the C 11 and O I lines with PDR models 
indicates a far-UV radiation field lO'^-lO' ^ times the mean 
interstellar field at the edge of Sgr B2 (|Goicoechea et al] 
L2OO4). This is consistent with the origin of water and OH in 
these photoactive layers. 

These lines are also particularly interesting towards 
Sgr B2 because the 63 /nm O i line and 158 lira C 11 line 
are seen in absorption in the line of sight clouds, and 
O I is self-absorbed at the velocity of Sgr B2 itself. This 
O I absorption was als o observed in the grating observations 
l|Baluteau et al.|[l997h , indicating a large fraction of the oxy- 
gen was in atomic form. The FP data confirm this, and show 
that not only is there absorption in the line of sight but the 
emission from Sgr B2 itself is completely cancelled by ab- 
sorpt ion in its outer envelope l|Lis et al.ll200ll : IVastel et al.l 
|2002| ). In contrast, emi ssion is observed throughout the rest 
of the extended cloud jBaluteau et al.ll 19971 : iLis et aLlliooH : 
IGoicoechea et al.l l2004'). 

The higher excitation O I 145 /im line is seen purely in 
emission at the velocity of Sgr B2. This is because the lower 
energy level of the transition is not sufficiently populated in 
the cool line of sight gas to show absorption. 

In the standard PDR view, there are 3 main layers 
going from atomic at the outside of the cloud, to where 
O I and C 11 coexist, to whe re all carbon is locked into CO 
IIHollenbach fc Tielens|[l999l '). In this last layer, O I coexists 
with CO until the O1/O2 transition which only occurs at 
much higher extinction levels. The high resolution FP ob- 
servations of O I and C II in absorption towards Sgr B2 are 
a very useful probe of the PDRs along its line of sight. The 
absorption lines are particularly useful as they give a direct 



measure of the column density in each cloud. IVastel et al.l 
(!2002') fitted the O I absorption using CO and H I observa- 
tions as a template to disentangle the different PDR layers 
that contain atomic oxygen in each cloud. First, H I absorp- 
tion observations were used to fit the component associated 
with the outer atomic layers and then the remaining O ab- 
sorption was associated with CO observations tracing the 
cold molecular cores. The results for each of the velocity 
features in the line of sight are shown in Table [51 

A similar method was used by iLis et ahl (|200ll ) , who 
compared the column density of O I in three velocity 
ranges to that of CO (after subtracting the O I compo- 
nent associated with H l) to give an apparent correlation of 
Ol/CO~9±1.3. No significant intercept was found for the 
relationship between Ol and CO, indicating the lack of an 
intermediate PDR layer where O I is present between the 
CO region and completely atomic skin traced by H I. This 
layer is predicted by PDR models to contain the transition 
CO/C i/Cii. 

However, the L03 observations of the C 11 line at 158 yLxn 
can be used to provide additional information, and using 
this additional lin e (as well as the Ol line at 145 /xm), 
IVastel et al.l (|2002l ) could separate the predicted excess of 
O I. They calculated a slightly lower ratio of Ol/CO in 
the cloud cores of 2.5±1.8, and an excess of O I between 
molecular and atomic regions of the clouds which indicated 
iV(Cl) = (2.4 ± 0.9) X 10" cm"^. This is in approximate 
agreeme nt with ground b ased observations of C I towards 
Sgr B2 l|Vastel et al.ll2002l ). These resuhs indicate ~70% of 
the oxygen is in atomic form and not locked into CO. Un- 
certainties in this analysis are large due to the spectral reso- 
lution of the LWS FP, and observations of C 11 and C I with 
Herschel HIFI will give a clearer picture of the different PDR 
layers. 

It has been proposed that the absorption of C 11 by fore- 
ground clouds as observed here could also explain observa- 
tions of bright galaxies which show a deficienc y in C 11 flux 
compared with their total FIR luminosity (jVastel et al.l 
|2002| ). The spectrum of the bright ultraluminous galaxy 
Arp220 shows a very similar spectrum to Sgr B2 with 
strong absorption by molecules (OH, H2O, CH, NH and 
NH3 ) as well as absorption by O I and weak emission by 
Cii l|Gonzalez- Alfonso et al.ll2004D . Gonzalez- Alfonso et al. 
found that they required some contribution from an absorb- 
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ing 'halo' to account for the observed O I absorption. How- 
ever, the deficiency in Cll/FIR could mostly be explained 
by the non-PDR component of the FIR continuum, although 
some effect from extinction in the halo probably also con- 
tributes. 

4-4-3 Ionised lines of oxygen and nitrogen 

These lines are associated with the warm ionised gas be- 
yond the edges of the PDRs, physically distinct from the 
region emitting/absor bing O I. The N ll line at 122 /jm 
was first observed by iRubin et a l. (19 89.) using the KAO 
towards G33 3 .6— 0. 2 (see also [Erickson et al.l l|l99ll ) and 
IColgan et all (|l993l )). N II is an important coolant in the 
ISM: the COBE satellite found that the two N ll fine struc- 
ture lines at 122 ^m and 205 /xm are the brightest in th e 
galaxy after the C ll line at 158 pim (Wright ct all Il99lf ). 
The N II emission is preferentially produced in lower density 
gas than the pair of O iii lines observed at 52 and 88 /xm. 

Low spectral resolution LWS grating rasters of the 
N II and O ill lines around Sgr B2 revealed a very extended 
component o f ionised gas, with aver age electron densities of 
~240 cm~^ (jGoicoechea et al]|2004l ). Detailed photoionisa- 
tion models of these lines show that the ionising radiation 
has an effective temperature of ^36,000 K and a Lym an con- 
tinuum flux of ^10^"" s"^ l|Goicoechea et all '2004V How- 
ever, the location and distribution of the ionising sources in 
the Sgr B2 envelope remains unclear, and will require much 
higher angular resolution observations. The LWS L04 mode 
was also used to o bserve N ii towards several other Galac- 
tic C entre sources (|Rodn'guez-Fernandez fc Martm-Pinta 
I2OO5I ). 

The N II line is particularly interesting as it has several 
hyperfine structure components that produce a noticeable 
broadening in the observed profile. This is due to the nu- 
clear spin of the ^'*N atom, resulting in 3 gr oups of com- 
pone nts which are separated by ~ 30 km s~^ (|Brown et al.l 
I1994I ). In order to fit the N 11 line, we fixed the relative con- 
tribution of each component to its expected line strength 
taken from Brown et al, (,1993), a^nd fixed the line width to 
51 km s~^ (which is the average of that observed for the 
other atomic lines above 70 /xm) and the centre velocity to 
be 65 km s^^. The strongest line component occurs at a 
wavelength of 121.8884 /im. Figure [T2l shows the result of 
the fit, giving a total N 11 fiux of 3.7 ± 0.7 x 10"^** W cm"^ 
(adopting a continuum level estimated from the LOl grating 
observation of 0.77 x 10" W cm ). Even accounting for 
the hyperfine components and a line width wider than the 
LWS FP spectral resolution at 122 (~34 km s~^), the fit- 
ted profile does not explain all the observed emission. There 
is some extra emission at a velocity of -1-200 km s~^ that is 
difficult to explain. A trace of this high velocity emission is 
also present in the O iii line at 51.8 pm (see Fig. Ilip . 

4.5 Global line properties and kinematics 

Figure [T3] shows how the results of the line fits vary across 
the survey range. In the upper plot showing the central ve- 
locity, there is a clear systematic error in the wavelength 
calibration for the non-prime data below 70 /im. This can 
be explained because the polynomial used to fit the relation- 
ship between the instrumental FP encoder setting and gap 
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Figure 12. The Nil line at 121.9 /^m. The solid line shows a fit 
including hyperfine splitting, with the relative strength of hyper- 
fine components set to their expected fine strengths (the thick 
bars show the location of the individual components). The line 
width at each component was fixed to 51 km s~^. 



between the FP plates for FPL was based on data above 
70 fj,m. This effect is also clear when the non-prime lines 
are directly compared with the equivale nt prime FPS obser- 
vation (e.g. the 53 OH lines - see iPolehampton et al.l 
l2005al ). A correction to the wavelength calibration has not 
been applied in the data reduction process for the survey, 
or in the results presented in Table [3] (however, it was taken 
into acc ount where data below 70 ixm were analysed in de- 
tail; e.g. IPolehampton et al.|[2005al ). 

In the wavelength range 70-196 /xm, the average veloc- 
ity over all lines is 62.6 km s~^ with an RMS of 5.5 km s~^. 
This was calculated by including only the absorption from 
Sgr B2 itself and ignoring the outlying lines above 80 km s~^ 
and below 50 km s~^. These outliers are: the H2O line at 
99.5 /xm (46 km s"^); the YifO line at 102.0 /xm (49 kms"^); 
the C3 /?(2) line at 152.3 /im (84 km s~^), already discussed 
in Sect. 14.3.31 the NH3 line at 165.7 ^im (82 km s"^); and the 
3 NH3 lines at 170 /xm which are blended, distorting the cen- 
tral velocity. Apart from these lines, the spread in velocity 
about the average is consistent with a study of the accuracy 
of the wavelength calibration using CO lines in the FP spec- 
trum of Orion - in that case, the line centres w ere always 
less t han 11 km s"'^ from the expected position (|Grv et al.l 
l2003h . 

We have also calculated the average velocity separately 
for each individual species. However, it is difficult to de- 
termine if differences are real or due to systematic ef- 
fects. In Fig. [13] there appears to be a trend with more 
lines above average >120 pim and more lines below aver- 
age <120 iJ,m. This is reflected in the average velocities 
because all the C3 lines are at wavelengths longer than 
120 /xm and most of the NH3 lines are in the shorter 
wavelength bracket. The average velocities for each species 
are, NH3 59.1±3.6 km s'^ NH2 63.5±0.6 km s"^ H2O 
60.8±4.2 kms"^ OH 60.2±8.4 kms'^ C3 66.7±4.8 km s"^ 
The atomic and ionic lines are generally higher in veloc- 
ity than the molecular lines, O ill (88 ^m) 72 km s~^; 
Ol (145 Mm) 75 km s"^; C 11 (158 ^m) 74 km s"\ 

Variation in the velocity of lines from different species 
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T able 6. Summary of the co lumn densities calculated for the line of sight features. The velocity ranges are based on t he H I me asurements 
of lGarwood fc DickevI lll989l). Only features <0 km s '^ were considered for the atomic lines. Columns 2-6 are from jVastel e t al. 20oJ), 
column 7 and 9 are from l lPolehampton et al.ll2005ah and column 8 is from IPolehampton et al"] ll2005lj) . Estimates of the uncertainties 
are given in those references. 
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was a lso found in the 330-355 GHz survey of ^Sut ton et al 



1991 ) and the 218-263 GHz survey of NummeUn et al 



200d ) . Sutton et al. found an average velocity of 60.6 km s 



for Sgr B2 M and Nummelin et al. found an average ve- 
locity of 61.6 km s^i, with smaller molecules in range 55- 
60 km s~^. However, at these wavelengths the lines are more 
likely to trace the material in the cores of the cloud rather 
than in the outer regions of the envelope as for our survey. 
The average lin e widths found in these mm surveys 
18.5 km s~^ JSutton et all Il99ll ) and 17 km s" 



iNummelin et al.ll2000t ). of order half the spectral resolution 
of our survey. The lower plot in Fig. 1 131 shows our measured 
line widths from Table O with the predicted instrumental 
profile width based on the measured resolving power (from 
Fig. lAip . The measured line widths cluster at or above this 
line, consistent with broadening caused by the width asso- 
ciated with the source. 

4.6 Isotopic species 

Several isotopic variants of the observed species were de- 
tected in the survey - these are due to the two species show- 
ing the strongest absorption in their main lines, H2O and 
OH. Observations of isotopologues can be used to determine 
intrinsic isotope ratios, if good estimates of the column den- 
sity of the main species can be made. However, if these ratios 
are already known, the isotopologue lines are also valuable 
in modelling as their optical depth is much lower than the 
main lines. 



4-6.1 Oxygen isotopes 

The ground state lines of both ^'^OH and ^*OH are observed 
in the survey, showing absorption in both Sgr B2 and the 
line of sight clouds. OH is particularly good for tracing the 
intrinsic oxygen isotopic ratio, because chemical fractiona- 
tion reactions that might act to distort their values (e.g. 
as occurs for ""^^C/^ ^C) are not thought to be important 
(|Langer et aL|[l984 ). In addition, optically thin lines from 
all three isotopologues were observed as part of the survey 
(for "'OH these were the cross ladder transitions). 

In order to calculate the i^O/^*0 ratio, the ground 
state line profile of i^OH w as compared with that of i*OH 
ijPolehampton et al.ll2005al l. The isotopic ratio in each ve- 
locity component of the lines was determined by carrying 



out a simultaneous fit of the i^OH 53 /.im and 79 /^m lines 
with the i*OH 120 /^m line. The fit was carried out in a 
similar way to that described for ^®0H in Sect. 14.3.21 us- 
ing H I observations as a basis for the line shapes of each 
velocity component. The results show values of i^0/^*0 
in the Galactic Disc that are broadly consistent with pre- 
vious measurements (in the range 360-540). However, in 
velocity components associated with the Galactic Centre, 
the ratio is higher than previous estimat es (3201,° com- 



pared with the standard value of 250; IWilson fc 
I1994D . The standard ratio for the Galactic Centre was de- 
rived fro m measurements of the A-doublet transitions of 
OH (e.g. IWiUiams fc Gardned Il98ll : I Whiteoak fc Gardner! 
Il98 j ). However, Buiarrabal et al.l l|l983l ) have suggested 



that these could be an underestimate due to excitation 
anomalies in the hyperfine levels caused by FIR rotational 
pumping. The survey results would appear to support this 
(the rotational lines at the resolution of the LWS FP are not 
affected by anomalies in the hyperfine levels). 

The only previous measurements of ^^OH in the ISM 
were m ade towards Sgr B2 via its m i crowave A-doublet tran - 
sitions l|Gardner fc WhiteoakI [l976l : iBuiarrabal et al.lll983D 
and the ISO observations are the first time its pure rota- 
tional transitions have been seen in space. Previous mea- 
surements of the value of i*0/^^0 in the Galactic Disk 
using CO indi cated a constant isotopic ratio of 3.65+0.15 
(Penzias| [l98ll ) and this agrees with the microwave obser- 
vations of ^^OH towards Sgr B2 (3.6+0.5; Buiarra bal et al.l 
[l_983). An alysis of the lines from th e survey has been pre- 
sented by IPolehampton et al.l l|2003l ). Due to the weakness 
of the ^^OH lines, they did not fit them directly, but used 
the ^*0H observations with a ratio of 3.5 (the weighted aver- 
age of all previous measurements in the Galactic Centre) to 
predict the line shape. This is complicated by the hyperfine 
splitting in ^'^OH (due to the spin of the ^^O nucleus) which 
causes a series of line components with similar separation 
to the velocity features along the line of sight. Using a ratio 
of i*OH/i^OH=3.5 shows good agreement at the velocity of 
Sgr B2 but some disagreement for the line of sight, appear- 
ing to underestimate the observed absorption (indicating a 
lower ratio is needed). However, this could also be due to 
uncertainties in the data r eduction distorting the line shape 
(IPolehampton et al.|[20o3 ). 

In comparison to OH, water is not a good tracer of 
the true oxygen isotopic ratio. This is because even though 
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Figure 13. Best fitting centre velocity (top) and line width (bot- 
tom) of identified features from Table [3] In the top plot the av- 
erage velocity in the range 70-196 fira is shown as the horizontal 
dashed line (62.7±5.6 km s~^). The vertical dashed line shows 
the edge of the nominal operating mode of FPL (70 ^m), be- 
low which the wavelength calibration has significant systematic 
errors. In the bottom plot the dashed line shows the predicted 
width based on the instrument resolving power. 



transitions of both H2'^0 and H2^0 were detected, they are 
highly optically thick and th e excitation of the ob served 
water line s is hiRhly complex |Cernicharo et al1l2006t ). The 
model of ICernicharo et af] (|2006l ) predicts that even the 
H2*0 line at 181.05 ^m (2i2-loi) is optically thick, with 
T ~ 5 (in addition, this line is blended with H3O"*'). The 
equivalent transition for H2^0 at 179.53 /^m has r ~ 10^- 
10^. The next strongest line is 22i-2io, and the model shows 
that r ~ 300 for H2^0. Therefore, water is much more dif- 
ficult to use as a probe of the isotopic ratio. 



4-6.2 Carbon isotopes 

The ratio ^^C/"C has been fo und to be 25±1 in Sgr B 2 de- 
rived from CO observations (iLanger fc PenziasI [l990l ) and 
so we might expect isotopologues of the carbon bearing 
species to be easily detected. However, ^ ^CH is not detected, 
with a lower limit on ^'^CH/^'^CH of 36 (|Polehampton et al.l 
l2005bl ). This is higher than the value derived from CO but 
could be due to isotopic fractionation reactions which would 



increase ^^CH/^^C H over the true ^^C/^^C isotopic ratio 
(|Langer et al.lll984l 'l. 

The other carbon bearing lines in the survey such as C3 
are too weak for their isotopologues to be detected. 

4.6.3 Nitrogen isotopes 

The isotopic ratio, ^"'N/^'^N, is about 400-500 in Sgr B2. 
This has been measured in radio observatio ns towards 
Sgr B 2 N in both emission and absorption lines (|Peng et al.l 
Il993h . Therefore, we would not expect to detect any lines 
from ^^NHs above the noise level in the survey. However, the 
ground state (Iq-Oo) transition of both '^^NHs and ^^NHs 
has been observed with the ODIN satellite towards Sgr B2 
marson et al.|[2005h . They observed strong absorption 
from ^■'NHs at all line of sight velocities and weak absorption 
from ^'^NHs at the velocity of Sgr B2. 



4.6.4 Deuterium 

The key deuterated ion, H2D"'", has been detected via its or- 
tho ground state rotatio nal transition 2i2-lii at 126.853 
ijCernicharo et al.|[2007h . This line appears blended with the 
NH2 221-I10 hne at 126.8014 /xm, but is clearly detected 
because the broad blended feature cannot be due to NH2 
alone. This is because all other detected NH2 lines arising 
from excited states are narrow and appear at the velocity 
of Sgr B2 itself. Removal of the contribution of NH2 shows 
that the H2 D"'" hne i s absorbed by all components along the 
line of sight (|Cernich aro ot al. 2007,). 

H2D"'" is an important molecule as it plays a crucial 
role in the deuter ation of other species in cold gas (e.g. 
iGerlich et al.|[2003) . The calculated column density at the 
velocity of Sgr B2 is 9x10^^ c m~'^ and in the line of sight 
clouds is (2.7-6.7) X 10" cm"^ l|Cernicharo et al.ll2007l ). 

The only other isotopic species containing deuterium 
that has possibly been detected in the survey is HD. A ~ 3(7 
emission feature is seen in co-added prime and non-prime 
data corresponding to the HD J=l-0 line at 112.0725 /xm 
( Polehampton ot al. 2002a). This leads to a D/H ratio in the 
Sgr B2 envelope in the range (0.2-11) x lO"'' . 

For other molecules, the formation process favours 
deuteration and so the abu ndances can be dramatic ally en- 
hanced over the D/H ratio dRoberts fc Mill ar 2000). 

HDO: The lowest transition in the range is 22i-2o2 at 
159.354 /xm, which is not detected. The equivalent transition 
to the saturated hne in H2O at 179 pim 2i2-loi occurs out- 
side the LWS range at 235 ^m. However, ot her transitions 
of HD O have been observed from the ground (jComito et al.l 
|2003| . and references therein) , and some of the higher energy 
transitions in the FIR have be en observed in the ISO spec- 
tral survey towards Orion KL l|Lerate et al.|[200^ '). 

OD: Chemical models predict that OD/OH should be 
enhanced over the D/H ratio. In di ffuse clouds OD/OH 
may be between 3x lO""* and 4x 10"^ (jCrosweU fc Dalgarnol 
Il985h . This would mean that the ground state rotational 
transition would be clearly detected towards Sgr B2. How- 
ever, only the cross ladder transitions from the ground state 
occur within the LWS range and these are too weak to be 
detected. The ground sta te transition occurs at 215.4 yLxa 
l|Brown fc Schubertlll98 j ). unfortunately falling in the un- 
observable gap between bands 5 and 6 of Herschel HIFI. 



20 E. T. Polehampton et al. 



NH2D: This has been detected towards th e Sgr B2 N 
hot c ore with a NH2D/NH3 ratio of 1.1x10"^ (|Peng et alj 
I1993I ). This is enhanced compared with the local D/H ratio 
but is an order of magnitude smaller than in hot core sources 
outside of the Galactic Centre. This ratio is too low to expect 
a detection in the survey. 

4.7 Unassigned features 

We have catalogued the unidentified features in the survey 
which have a strength greater than ~3cr above the RMS 
noise in data. All are seen in absorption and there are no 
unassigned emission lines (> 3a"). Possible assignments for 
these absorption features were checke d using the Cologn e 
Database for Molecular Spectr oscopy |Miiller et al.ll2005h . 
the JPL spectral line database (iPickett et al.lll998l ) and the 
catalogue of molecular transitions maintained by J. Cer- 
nicharo. Table [7] gives a list of these features with their 
fitted line widths and depths, which were determined in a 
similar way to the assigned lines: the line centre, width and 
depth were free parameters in a Lorentzian fit. The wave- 
lengths were calculated assuming that the lines were centred 
at 65 km s~^. All the wavelengths for unidentified features 
mentioned in the following text have been corrected by this 
velocity. In addition to the strong features listed in Table [3 
there are many more possible lines with a significance < Scr. 

The features in Table [7] were fitted in the binned data, 
but have also been checked in the raw unaveraged data to 
make sure they were not due to glitches, thin patches in 
the data where glitches had been removed or bad matching 
where mini-scans overlap. Such spurious features are easy to 
spot in the unaveraged data if they are examined mini-scan 
by mini-scan. All the features in Table [7] are considered to 
be real, and in most cases appear in at least 2 independent 
observations. The level of each detection above the RMS 
noise in the adjacent continuum is also shown. The line pro- 
file of each unidentified line are shown in detail in Fig. 1151 
Some of the features are narrow, with a width close to the 
FP resolution, and some are broad indicating possible con- 
tribution from the line of sight clouds (or multiple blended 
line components). The 139.12 /im feature appears to have 2 
distinct components. 

The feature at 180.338 ^m is close to the 2i2-loi line 
of H2^0 at 180.3302 /^m, but if this is really the true iden- 
tification, the line would peak at 79 km s~^ rather than the 
expected 65 km s~^. 

There are three features in the list that are particu- 
larly notable because they clearly appear in multiple obser- 
vations, including ones made on different detectors: these 
are at 153.47 /^m, 175.07 /im and 184.27 /^m. The feature 
at 175.07 /im is particularly strong (~6(t), clearly showing 
a broad profile, similar in shape to the lines of CH and OH 
due to absorption along the whole line of sight to Sgr B2. 
Figure [HI shows the line shape compared with that of OH. 
This indicates that the line must be due to a low energy 
transition, probably from the ground state. In all the other 
molecules observed there are no higher energy transitions 
that show line of sight absorption. However, we have not 
been able to identify any candidates that fit this descrip- 
tion. The closest matches found in the line catalogues are 
all too high in energy: HsO^ 6g'-6^ at 175.0633 ^m; H2O 
5o5-432 at 174.9200 ^.m. 



One possible group of species that could be responsible 
for some of the unidentified features is the ions of the neu- 
tral hydride species that were observed in the survey. These 
species may occur in the external layers of the Sgr B2 en- 
velope where the volume density is low (in the inner parts 
they would be destroyed by reaction with H2). In the follow- 
ing paragraphs we present examples of some of the hydride 
cations that could be thought to be detected. However, none 
of these molecules match with any of our unassigned lines, 
showing that further study and lab work will be needed to 
identify the true carriers. 

Rotational transitions of the CH^ ion have been de- 
tected in emission in t he planetary nebula NGC 7027 
(|Cernicharo et al]|l997bl ). However, its lowest energy tran- 
sition in our survey range, J = 2-1, at 179.611 pm overlaps 
with the much broader saturated absorption line due to wa- 
ter and so cannot be separated. 

The OH+ ion is isoelectronic to the NH radical and 
could be p roduced in the outer layers of interstellar clou ds or 
in shocks (|de Almeida fc Singhlll98ll : Ide Almeidalll990l . and 
references therein). The fundamental transition, A'^ = 1-0, 
occurs at 329.8 /xm U = 0-1), 308.5 A ^m (J = 2-1) and 
290.2 /xm (J = 1-1) (jSekoov et al.lll985l ') outside the range 
of our survey. However, there is also a cluster of features 
that could match with the broad U-line at 153.47 /xm, due 
to the = 2-1, J = 1-2 transition (Cologne Database 
for Molecular Spectroscopy). Unfortunately, this transition 
has a very low line strength, and it would be much more 
likely to detect the stronger feature due to the N = 2-1, 
J = 3-2 transition at 152.9897 ^m. This is not detected, 
ruling out any identification of OH""" at 153.47 /im. A much 
better chance of detecting OH"^ will be provided by the HIFI 
instrument on board the Herschel satellite, as it will be able 
to observe the ground state lines. 

The NH"*" io n is isoelectronic to CH and 
Ide Almeida fc Singhl (|l982l ) have proposed that pho- 
toionisation, charge-transfer, radiative association and 
ion-molecule reactions in the gas phase may lead to a 
significant concentration of NH"*" in the ISM. Its ground 
state lines {N = 1-1, J = 3/2-1/2) oc cur outside our 
surve y range at 294.14 /xm and 296.08 /xm (|Verhoeve et al.l 
Il986l l. but the second rotational transitions, A'^ — 2-1, 
J = 5/2-3/ 2, do occur within the ran^ e at 160.01 /im and 
159.70 /im (|Kawaguchi fc Amand ll988l. However, these do 
not correspond to any of the unassigned lines. 

Calculations for CH^ show that its low energy ro- 
tational transitions all occur within the survey range. 
However, accurate measurements of its rotational levels 
have not been made in the laboratory. The best mea- 
surements so far have be en made by photoelectron spec- 
troscopy (Wi llitsch fc Me rkt 2003) and by near-IR spec- 
troscopy (jCottfried fc Okal |200J). Using the parameters 
from these measurements, we calculate that the ground state 
lii-Ooo transitions should be at ~128.9 /xm (J = 3/2-1/2) 
and ~ 134.6 /xm ( J = 1/2-1/2). The next higher energy tran- 
sitions are 2ii-loi at ^109. 6 /im and lio-loi at ~157.2 /im 
and 164.5 fim. There is considerable uncertainty in these pre- 
dictions (up to 2 /im) , in particular because the ^-rotational 
constant and the spin-rotation splittings are very poorly 
determined or not even known, and more lab experiments 
are needed before a firm identification could be made. The 
molecule is predicted to have a reasonably sized dipole mo- 
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Figure 14. The line shape of the unidentified feature at 
175.07 /im (above) compared with the ^'^1/2-^^3/2 J=l/2-3/2 
fine of OH at 79 lira (below). The OH line has been scaled to 
match the unidentified feature. 



ment, calcul ated to be 0.629 D by ([psmann et aLlll997h and 
0.701 D by (|Brinkmann et al.ll200a V 

The H2O''' ion has been observed in optic al emission 
spectra of comets (e.g. ICremonese et all |2007^ . However, 
its ground state rotational transitions, In- Opo, occur at 
268.84 ^im and 263.08 ^im (|Murtz et al.lll998l '). outside our 
survey range. The next higher energy transitions are 221-I10 
at 105.74 ^m and 104.72 /xm and 220-I11 at 100.91 ^m and 
99.88 ^m l|Mflrtz et al.lll998l i. None of these are detected. 

In addition to the lines from simple molecules, there 
could be featur es from the bendin g modes of carbon clus- 
ters and PAHs (|Mulas et al.l 120061 ). In emission, these lines 
show a band structure due to the Q-branch lines (which 
would be resolved by the LWS FP) and a broader 'grass'. 
No emission features matching the predicted band structure 
of PAH molecules has been observed in the survey, however, 
it is possible that the broader component could exist and 
was removed in the reduction process. This is because very 
broad features in the continuum in the FP data are impos- 
sible to distinguish from artifacts of the reduction process 
(see Sect. 13.4]) and so were removed. However, in absorption 
by cold gas, only the ground state Q-branch line would be 
observed and the profile could be narrow. In order for this 
absorption to be detected, a molecule with high transition 
intensity would be required. Concrete identification of a U- 
line with a single PAH molecule would be very difficult in 
this case, and require further laboratory measurements of 
the candidate species. 



5 SUMMARY 

We have presented the data from a complete spectral sur- 
vey of Sgr B2 between 47 and 196 /im, observed with the 
ISO LWS FP spectrometer. The spectral resolution was 30- 
40 km s~^ across the range and the signal-to- noise achieved 
in the line to continuum spectrum was ~100. A total of 95 



Table 7. Unidentified Lines in the survey data. The wavelengths 
given were calculated assuming the line centres are at 65 km . 
Approximate widths and line depths were derived by fitting the 
lines with a Lorentzian profile in the same way as the assigned 
lines in Table [3] The 4th column gives the detection above the 
RMS noise in the surrounding continuum. The lines in bold are 
the most secure (i.e. detected in multiple observations). 
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0.836 
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184.274 


60 
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lines have been identified in the spectrum, with a further 11 
unassigned features. The main results can be summarised as 
follows; 

• The survey is dominated by molecular lines in absorp- 
tion due to NH3, NH2, NH, H2O, OH, H3O+, CH, CH2, C3, 
HF and H2D^ (and in emission for OH at 163 /im). The 
ground state lines of OH, CH, CH2, ortho-H20, H3O''' and 
H2D^ show absorption due to the entire line of sight as well 
as Sgr B2 itself. Isotopic lines are observed from ^®OH, ^'^OH 
and }lfO. 

• Atomic and ionic lines of O I, O ill, C 11, N 11 and 
N III are seen in emission at the velocity of Sgr B2 itself, 
with the exception of the O I 63 /xm line which is completely 
self-absorbed. This line and the C 11 158 /im line show ab- 
sorption due to the line of sight clouds. 

• The species with greatest column density in the survey 
is H2O, with iV(Sgr B2) ~ 9 X 10^^ cm'^ 

• The average velocity for all the assigned lines which only 
show absorption due to Sgr B2 is 63±6 km s~^. There ap- 
pears to be some variation in the velocity of different species 
but it is difficult to separate this from systematic effects. 

• Several unassigned features remain, all of which are in 
absorption. These could be due to simple molecules or bend- 
ing modes of carbon chains/clusters. 

This survey provides a basis for future follow up at 
higher spectral and spatial resolution (e.g. with HIFI, or 
instruments on SOFIA). The planned Herschel guaranteed 
time spectral survey of Sgr B2 with HIFI and PACS will 
overlap with this survey above 57 /^m, and extend it into 
the sub-mm. It will be interesting to check the unassigned 
features in the overlap region, particularly the line at 175 /xm 
which shows a very similar profile to the ground state lines 
of OH and CH. 



ACKNOWLEDGEMENTS 

We would like to thank our referee for useful suggestions, 
and Christine Joblin, Holger Miiller and Mercedes Lerate 



22 E. T. Polehampton et al. 




-300 -200 -100 100 200 300 -300 -200 -100 100 200 300 -300 -200 -100 100 200 300 
"lsr (km s"') i;lsr (km s"') ij^s„ (km s'') 

Figure 15. Unidentified lines in the survey data plotted versus velocity (assuming they are centred at 65 km s~^). Note that the y-axis 
has a different scale for each of the three plots. 

for helpful discussions. JRG was supported by an individ- 
ual Marie Curie fellowship, contract MEIF-CT-2005-515340. 

The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and 
Data Centres. Contributing institutes are CESR, IAS, IPAC, 
MPE, RAL and SRON. LIA is a joint development of the 
ISO-UWS Instrument Team at the Rutherford Appleton 
Laboratory (RAL, UK - the PI Institute) and the Infrared 
Processing and Analysis Center (IPAC/Caltech, USA). 



REFERENCES 

Adamkovics, M., Blake, G. A., & McCall, B. J. 2003, ApJ, 
595, 235 

Baluteau, J.-P., Cox, P., Cernicharo, J., et al. 1997, A&A, 
322, L33 

Bekooy, J. P., Verhoeve, P., Meerts, W. L., & Dymanus, A. 

1985, J. Chem. Phys., 82, 3868 
Brinkmann, N. R., Richardson, N. A., Wesolowski, S. S., 

Yamaguchi, Y., & Schaefer III, H. F. 2002, Chem. Phys. 

Lett., 352, 505 

Brown, J. M. & Schubert, J. E. 1982, J. Molec. Spectrosc, 
95, 194 

Brown, J. M., Schubert, J. E., Evenson, K. M., & Radford, 

H. E. 1982, ApJ, 258, 899 
Brown, J. M., Varberg, T. D., Evenson, K. M., & Cooksy, 

A. L. 1994, ApJ, 428, L37 
Bujarrabal, V., Cernicharo, J., & Guelin, M. 1983, A&A, 

128, 355 

Ceccarelh, C, Baluteau, J.-P., Walmslcy, M., et al. 2002, 

A&A, 383, 603 
Cernicharo, J., Goicoechea, J. R., & Caux, E. 2000, ApJ, 

534, L199 

Cernicharo, J., Goicoechea, J. R., Pardo, J. R., & Asensio- 

Ramos, A. 2006, ApJ, 642, 940 
Cernicharo, J., Lim, T., Cox, P., et al. 1997a, A&A, 323, 

L25 



Cernicharo, J., Liu, X.-W., Gonzalez-Alfonso, E., et al. 

1997b, ApJ, 483, L65 
Cernicharo, J., Polehampton, E., & Goicoechea, J. R. 2007, 
ApJ, 657, L21 

Clegg, P. E., Ade, P. A. R., Armand, C, et al. 1996, A&A, 

315, L38 

Colgan, S. W. J., Haas, M. R., Erickson, E. F., et al. 1993, 
ApJ, 413, 237 

Comito, C, Schilke, P., Gerin, M., et al. 2003, A&A, 402, 
635 

Cooksy, A. L., Blake, G. A., & Saykally, R. J. 1986, ApJ, 

305, L89 

Cox, P., Baluteau, J.-P., Roelfsema, P. R., & Swinyard, B. 

1999, in Astrophysics with Infrared Surveys: A Prelude 

to SIRTF, ASP Conference Series 177, 275 
Cremonese, G., Capria, M. T., & De Sanctis, M. C. 2007, 

A&A, 461, 789 
Croswell, K. & Dalgarno, A. 1985, ApJ, 289, 618 
Cummins, S. E., Linke, R. A., & Thaddeus, P. 1986, ApJS, 

60, 819 

Davidson, S. A., Evenson, K. M., & Brown, J. M. 2001, 

ApJ, 546, 330 

Davis, G. R., Furniss, I., Towlson, W. A., et al. 1995, Appl. 
Optics, 34, 92 

de Almeida, A. A. 1990, Rev. Mexicana Astron. Astrof., 
21, 499 

de Almeida, A. A. & Singh, P. D. 1981, ApJ, 95, 383 
de Almeida, A. A. & Singh, P. D. 1982, A&A, 113, 199 
Emery, R. J., King, K. J., Swinyard, B. M., & Church, S. E. 

1993, SPIE, 1946, 272 
Erickson, E. F., Colgan, S. W. J., Simpson, J. P., et al. 

1991, ApJ, 370, L69 
Evenson, K. M., Jennings, D. A., Brown, J. M., et al. 1988, 

ApJ, 330, L135 
Flower, D. R., Pineau des Forets, G., & Walmsley, C. M. 

1995, A&A, 294, 815 
Priedel, D. N., Snyder, L. E., Turner, B. E., & Remijan, A. 

2004, ApJ, 600, 234 



The ISO LWS spectral survey towards Sgr B2 23 



Gardner, F. F. & Whitooak, J. B. 1976, MNRAS, 176, 57P 
Garwood, R. W. & Dickey, J. M. 1989, ApJ, 338, 841 
Gerlich, D., Herbst, E., & Roueff, E. 2002, P&SS, 50, 1275 
Giesen, T. F., Van Orden, A. O., Cruzan, J. D., ct al. 2001, 

ApJ, 551, L181 
Goicoechea, J. R. & Cernicharo, J. 2001, ApJ, 554, L213 
Goicoechea, J. R. & Cernicharo, J. 2002, ApJ, 576, L77 
Goicoechea, J. R., Rodriguez- Fernandez, N. J., & Cer- 
nicharo, J. 2004, ApJ, 600, 214 
Goldsmith, R F., Lis, D. C, Hills, R., & Lasenby, J. 1990, 
ApJ, 350, 186 

Goldsmith, P. F., Lis, D. C, Lester, D. F., & Harvey, P. M. 

1992, ApJ, 389, 338 
Gonzalez- Alfonso, E., Smith, H. A., Fischer, J., & Cer- 
nicharo, J. 2004, ApJ, 613, 247 
Gottfried, J. L. & Oka, T. 2004, J. Chem. Phys., 121, 11527 
Greaves, J. S. & Nyman, L. 1996, A&A, 305, 950 
Greaves, J. S. & Williams, P. G. 1994, A&A, 290, 259 
Griffin, M. J. & Orton, G. S. 1993, Icarus, 105, 537 
Cry, C, Swinyard, B., Harwood, A., ct al. 2003, ISO Hand- 
book Volume III (LWS), Version 2.1, ESA SAI-99-077/Dc 
Hinkle, K. W., Keady, J. J., & Bernath, P. F. 1988, Science, 
241, 1319 

Hjalmarson, A., Bergman, P., Biver, N., et al. 2005, Ad- 
vances in Space Reseaxch, 36, 1031 

HoUenbach, D. J. & Tielens, A. G. G. M. 1999, Reviews of 
Modern Physics, 71, 173 

HoUis, J. M., Jewell, P. R., & Lovas, F. J. 1995, ApJ, 438, 
259 

Hiittemeister, S., Wilson, T. L., Henkel, C, & Mauers- 

bcrgcr, R. 1993, A&A, 276, 445 
Johns, J. W. C. 1985, J. Opt. Soc. Am. B, 2, 1340 
Kawaguchi, K. & Amano, T. 1988, J. Chem. Phys., 88, 

4584 

Langer, W. D., Graedel, T. E., Frerking, M. A., & Armen- 

trout, P. B. 1984, ApJ, 277, 581 
Langer, W. D. & Pcnzias, A. A. 1990, ApJ, 357, 477 
Leeks, S. J. 2000, PhD thesis. Queen Mary and Westfield 

College, University of London 
Lerate, M. R., Barlow, M. J., Swinyard, B. M., ct al. 2006, 

MNRAS, 370, 597 
Lim, T. L., Hutchinson, G., Sidher, S. D., et al. 2002, SPIE, 

4847, 435 

Lim, T. L. & Swinyard, B. M. 2002, in The Cahbration 
Legacy of the ISO Mission, ESA SP-481, 103 

Lis, D. C. & Goldsmith, P. F. 1990, ApJ, 356, 195 
Lis, D. C, Keene, J., Phillips, T. G., et al. 2001, ApJ, 561, 
823 

Lloyd, C. 2002a, in The Cahbration Legacy of the ISO 

Mission, ESA SP-481, 403 
Lloyd, C. 2002b, in The Cahbration Legacy of the ISO 

Mission, ESA SP-481, 399 
Lloyd, C, Lerate, M. R., & Grundy, T. W. 2003, The 

LWS LOl Pipeline, version 1, available from the ISO Data 

Archive at http://www.iso.vilspa.esa.es/ida/ 
Lugtcn, J. B., Stacey, G. J., & Genzel, R. 1986, BAAS, 18, 

1007 

Lyu, C.-H., Smith, A. M., & Bruhweiler, F. C. 2001, ApJ, 

560, 865 

Martin-Pintado, J., de Vicente, P., Wilson, T. L., & John- 
ston, K. L. 1990, A&A, 236, 193 
Matsushima, F., Nagase, H., Nakauchi, T., Odashima, H., 



& Takagi, K. 1999, J. Mol. Spcctrosc, 193, 217 
Matsushima, F., Odashima, H., Iwasaki, T., Tsunekawa, 

S., & Takagi, K. 1995, J. Mol. Struct., 352/353, 371 
Millar, T. J., Bennett, A., Rawlings, J. M. C, Brown, P. D., 

& Charnlcy, S. B. 1991, A&AS, 87, 585 
Morino, I. & Kawaguchi, K. 1997, J. Molec. Spectrosc, 182, 

428 

Morino, I., Odashima, H., Matsushima, F., Tsunekawa, S., 

& Takagi, K. 1995, ApJ, 442, 907 
Mulas, G., Malloci, G., Joblin, C, & Toublanc, D. 2006, 

A&A, 456, 161 
Miiller, H. S. P., Schloder, F., Stutzki, J., & Winnewisser, 

G. 2005, J. Mol. Struct., 742, 215 
Miirtz, P., Zink, L. R., Evcnson, K. M., & Brown, J. M. 

1998, J. Chem. Phys., 109, 9744 
Neufeld, D. A., Ashby, M. L. N., Bergin, E. A., et al. 2000, 

ApJ, 539, Llll 
Neufeld, D. A., Bergin, E. A., Melnick, G. J., & Goldsmith, 

P. 2003, ApJ, 590, 882 
Neufeld, D. A., Kaufman, M. J., Goldsmith, P. F., HoUen- 
bach, D. J., & Plume, R. 2002, ApJ, 580, 278 
Neufeld, D. A., Zmuidzinas, J., Schilkc, P., & Phillips, T. G. 

1997, ApJ, 488, L141 
Nolt, I. G., Radostitz, J. V., DiLonardo, G., et al. 1987, J. 

Molec. Spectrosc, 125, 274 
Nummelin, A., Bergman, P., Hjalmarson, A., et al. 1998, 

ApJS, 117, 427 
Nummelin, A., Bergman, P., Hjalmarson, A., et al. 2000, 

ApJS, 128, 213 
Osmann, G., Bunker, P. R., Jensen, P., & Kraemer, W. P. 

1997. Chem. Phys., 225, 33 
Peng, Y., Vogcl, S. N., & Carlstrom, J. E. 1993, ApJ, 418, 

255 

Penzias, A. A. 1981, ApJ, 249, 518 

Pettersson, S.-G. 1982, Physica Scripta, 26, 296 

Pickett, H. M., Poyntcr, R. L., Cohen, E. A., et al. 1998, 

J. Quant. Spectrosc. & Rad. Transfer, 60, 883 
Plume, R., Kaufman, M. J., Neufeld, D. A., et al. 2004, 

ApJ, 605, 247 

Polehampton, E. T. 2002, PhD thesis, Oxford University 

Polehampton, E. T., Baluteau, J. -P., Ceccarelli, C, Swin- 
yard, B. M., & Caux, E. 2002a, A&A, 388, L44 

Polehampton, E. T., Baluteau, J.-P., & Swinyard, B. M. 
2005a, A&A, 437, 957 

Polehampton, E. T., Brown, J. M., Swinyard, B. M., & 
Baluteau, J.-P. 2003, A&A, 406, L47 

Polehampton, E. T., Menten, K. M., Briinken, S., Win- 
newisser, G., & Baluteau, J.-P. 2005b, A&A, 431, 203 

Polehampton, E. T., Swinyard, B. M., Sidher, S. D., & 
Baluteau, J.-P. 2002b, in The Calibration Legacy of the 
ISO Mission, ESA SP-481, 435 

Roberts, H. & Millar, T. J. 2000, A&A, 361, 388 

Rodriguez-Fernandez, N. J. & Martm-Pintado, J. 2005, 
A&A, 429, 923 

Rubin, R. H., Haas, M. R., Erickson, E. F., Simpson, J. P., 
& Colgan, S. W. J. 1989, BAAS, 21, 1213 

Rudy, D. J., Muhleman, D. O., Bcrgc, G. L., Jakosky, 
B. M., & Christensen, P. R. 1987, Icarus, 71, 159 

Schmuttenmaer, C. A., Cohen, R. C, Pugliano, N., et al. 
1990, Science, 249, 897 

Sidher, S. D., Griffin, M. J., Davis, G. R., et al. 2000, Icarus, 
147, 35 



24 E. T. Polehampton et al. 



Spitzer, L. 1978, Physics of the IntersteUar Medium (New 

York: John Wiley and Sons) 
Stacey, G. J., Lugten, J. B., & Genzel, R. 1987, ApJ, 313, 

859 

Sternberg, A. & Dalgarno, A. 1995, ApJS, 99, 565 
Storey, J. W. V., Watson, D. M., & Townes, C. H. 1981, 
ApJ, 244, L27 

Sturm, E., Bauer, O. H., Lutz, D., et al. 1998, in Astro- 
nomical Data Analysis Software and Systems VII, ASP 
Conference Series 145, 161 

Sutton, E. C, Jaminet, P. A., Danchi, W. C., & Blake, 
G. A. 1991, ApJS, 77, 255 

Swinyard, B. M., Burgdorf, M. J., Clegg, P. E., et al. 1998, 
SPIE, 3354, 888 

Turner, B. E. 1991, ApJS, 76, 617 

van der Tak, F. F. S., Belloche, A., Schilke, P., et al. 2006, 

A&A, 454, L99 
van Dishoeck, E. F., Jansen, D. J., Schilke, P., & Phillips, 

T. G. 1993, ApJ, 416, L83 
Varberg, T. D. & Evenson, K. M. 1993, J. Molec. Spec- 

trosc, 157, 55 

Vastel, C., Polehampton, E. T., Baluteau, J.-P., et al. 2002, 
ApJ, 581, 315 

Vejby-Christensen, L., Andersen, L. H., Heber, O., et al. 

1997, ApJ, 483, 531 
Verhoeve, P., Ter Meulen, J. J., Meerts, W. L., & Dymanus, 

A. 1986, Chem. Phys. Lett., 132, 213 
Viti, S., Roueff, E., Hartquist, T. W., Pineau des Forets, 

G., & Williams, D. A. 2001, A&A, 370, 557 
Whiteoak, J. B. & Gardner, F., F. 1981, MNRAS, 197, P39 
Williams, D. R. & Gardner, F. F. 1981, PASP, 93, 82 
Willitsch, S. & Merkt, F. 2003, J. Chem. Phys., 118, 2235 
Wilson, T. L., C, H., & Hiittemeister, S. 2006, A&A, 460, 

533 

Wilson, T. L. & Rood, R. T. 1994, ARA&A, 32, 191 
Wright, E. L., Mather, J. C, Bennett, C. L., et al. 1991, 
ApJ, 381, 200 

Zink, L. R., Evenson, K. M., Matsushima, F., Neils, T., & 

Robinson, R. L. 1991, ApJ, 371, L85 
Zmuidzinas, J., Blake, G. A., Carlstrom, J., Keene, J., & 

Miller, D. 1995, ApJ, 447, L125 



APPENDIX A: DETAILED DATA REDUCTION 
Al Dark current and stray light 

Careful determination of the detector dark signals was very 
important for the accurate calibration of relative line depths 
in the survey. It is a particular problem for Sgr B2 due to 
the strength of the source and the fact that it was observed 
slightly off-axis. Scattering of off-axis rays (due to the par- 
tiall y reflecting sub strate of one of the LWS mirrors, M2; 



Table Al. Dark current and stray light values for FPS and FPL. 
The nominal dark currents are given for comparison - these are the 
combination of results made throughout the mission us ing three 
independent means of measurement given by 'Gry et al] pOOij ) . 



Grv et al.|[2003 l caused a series of reflections within the 



body of the instrument. These reached the detectors outside 
of the main beam without passing through the grating or 
FP and hence produced a background signal independent of 
wavelength. This was particularly important during FP ob- 
servations because the radiation not transmitted by the FP 
(> 99% of the incident flux) was reflected back towards M2. 
Stray light within the beam (wavelength dependant) was not 



Detector 


Nominal dark 


FPS 


FPL 




(10-1^ A) 


dark-f stray light 
(10-1*' A) 


dark+stray light 
(10-l» A) 


SWl 


4.96±0.54 


(a) 




SW2 


2.08±0.43 


(a) 


3.0('') 


SW3 


2.20±0.21 


2.48±0.09 


3.3±0.25(*) 


SW4 


1.18±0.34 


1.33±0.04 


1.6±0.1 


SW5 


1.56±0.24 


1.80±0.03 


1.9±0.2 


LWl 


2.50±0.29 


2.8±0.1 


2.7±0.l('=) 


LW2 


0.07±0.27 


0.60±0.07 


0.46±0.03 


LWS 


0.53±0.39 


2.2±0.4 


1.41±0.2 


LW4 


1.76±0.42 


4.1±0.6 


2.52±0.02 


LWS 


1.21±0.25 


1.8±0.1 


1.43±0.02 



Nominal value assumed 

Estimated from comparison of lines measured with both FPs 
^'^^ Sp ecial modelling performed for LWl - see iPolehamptonI 
ll2002f) . 



important because on-axis scattering was minimised by the 
design of the instrument. 

In order to characterise the dark + stray light towards 
Sgr B2, we used data recorded on non-prime detectors where 
the separation between FP orders was enough that no or- 
der occurred within the grating response function. At wave- 
lengths where the FP order separation was sufficiently large 
that any contribution from the FP profile wings was negligi- 
ble, the source radiation in the main beam was completely 
blocked. In this way, the dark -I- stray light signal could be 
accurately determined for most detectors. For the remaining 
detectors, it was estimated where possible by carefully com- 
paring lines that had been measured using both FPs. For 
the detector LWl using FPL, special modelling was used 
to estimate the dar k signal and this is described in full by 
IPolehamptoiil \2m% . The final values adopted are shown in 
Table Ell 

A similar method has been adopted to calculate the 
dark -I- stray light signal for th e ISO LWS FP survey to- 
wards Orion (jLerate et al.|[200^ ). A comparison of the re- 
sults shows that the contribution from stray light is much 
higher for Orion. This probably reflects both the differences 
in source morphology in the beam, and the different spectral 
shapes - Orion is very bright at shorter wavelengths with a 
peak in thermal emission <47/im. 



A2 Characterising the instrumental response and 
throughput 

The relative spectral response function (RSRF) of each 
LWS detector was determined in the grating mode us- 
ing observations of Uranus, and monitored throughout the 
mission using this and several other sources. Uranus was 
used as the primary calibrator because it is bright, point- 
like in the LWS beam and almost spectrally featureless 
and has a well establi shed atmospheric emission model (e.g 
iGrifHn fc Ortonlll993l ). The RSRFs were obtained by remov- 
ing a model of the Uranus spectrum from the observations. 
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Figure Al. Top: Resolving power of FPS and FPL (defined 
as A/AAfwhm) from ground- based and in-fiight measurements. 
Bottom: Third order polynomial fit to the total measured FP 
throughput for FPS and FPL (defined as the product of FP trans- 
mission efficiency and resolution element equivalent width). 



giving the wavelength dep e ndant response for eac h detec- 
tor band (|Grv et all booj : iLim fc SwinvardI l2002i '). These 
are dominated by the transmission of the band-pass filters 
placed in front of each detector. An additional correction 
was then applied to reference the absolute responsivity of 
each detector to its value when the RSRF calibration ob- 
servations were made. This was determined using standard 
illuminator flashes before and after each observation. 

The spectral response function of the LWS grating is 
also important for the calibration of FP data as it was used 
as an order selector for each FP and its shape is imprinted on 
all FP mini-scans. We used this fact to recover the grating re- 
sponse from non-prime FP data measured during the survey. 
This is possible because the FP scanned different regions of 
the profile in different non-prime observations. These were 
then recombined and the profile reconstructed. Data con- 
taining strong lines were excluded from the analysis and the 
Sgr B2 continuum spectrum was removed using the LOl grat- 
ing observation. The resulting data were fitted using a the- 
oretical profile (the convolution of a Gaussian profile for the 
grating with the FP Airy function) with its width as a free 
parameter. This gave a good fit and sets of grating profiles 
were built up for different grating angle ranges. This accu- 
rate knowledge of the grating response in FP mode allowed 



a further 'side order correction' to be made in the calibra- 
tion, accounting for neighbouring FP orders that may have 
been trans mitted th rough the wings of the grating response 
profile (see lGrv et~al. 2003: Polehampton ct al. 2002b). 

The remaining calibration needed for FP data is due 
to the properties of the FPs themselves. During the ground 
testing of the LWS, the transmission efficiency and spec- 
tral response of each FP were measur ed at a number of 
spot wavelengths using a FIR laser (jEmerv et al ] Il993l : 
iDavis et al.l 1 19951 ). However, in flight, it was difficult to 
confirm these two quantities independently of each other 
due to the lack of sufficiently well characterised bright nar- 
row lines in astronomical sources. Only a few such sources 
(G36.3-H0.7, IC3568, NGC7023, NGC7027, GO.6-0.6) were 
observed through the mission and the resulting values for the 
FP spectral resolving power versus wavelength are shown 
in the top panel of Fig. lAll (note that the LWS handbook 
only contains measurements of re solving power ma de during 
ground testing of the instrument; iGrv et aLll2003l ). Detailed 
fitting of the lines showed that the FP response function 
could be accurately described by the theoretical Airy profile 
with coefficient of finesse set by the resolving power. In the 
case of NGC7023, the lines were sufficiently narrow to allow 
a direct comparison with the grating leading to a determi- 
nation of the FP transmission at two wavelengths (8.3% at 
63.2 for FPS and 3.5% at 157.7 ^.m for FPL). 

Due to the limited wavelength coverage of these flight 
measurements, only the product of transmission efficiency 
and spectral resolution element was used in the standard cal- 
ibration procedure. This quantity, the FP 'throughput', was 
easier to measure in practice by observing a bright source 
with known continuum spectrum by scanning the grating 
whilst keeping the FP gap fixed. This produces a convolution 
of the FP Airy function with the grating response, the peaks 
of which trace the FP throughput. These observations were 
carried out towards Mars and the source continuum removed 
using the Martian model develope d by Rudv et akl (|l987l ) 
adapted for the LWS spectral range IjSidher et al ]|2000l) . The 
proce dure has been de scribed previously in the LWS hand- 
book (|Grv et all [20031 ) but was improved and extended to 
wavelengths outside the nominal range of each FP to allow 
non-prime data from the spectral survey to be calibrated. In 
addition, an equivalent dataset was observed towards Sgr B2 
which shows reasonable agreement with the Sgr B2 LOl spec- 
trum when processed wi th the extended throughpu t correc- 
tion derived from Mars (jPolehampton et al. l l2002bD . The fi- 
nal throughput results are shown in the bottom panel of 
Fig. lAll The final transmission efficiency of each FP can 
be determined by combining both plots in Fig. lAll and this 
gives excellent agreement with the two points measured to- 
wards NGC7023. 

Figure IXT1 shows that there is an apparent break in the 
throughput where the coverage by detector LWl ends and 
LW2 begins. One possible explanation for this discontinuity 
is that detectors SWl-LWl were mounted on a physically 
separate structure to LW2-LW5 and may have had different 
optical alignment with respect to the pupil. Another possi- 
bility is due to the detectors themselves - LW2-LW5 all used 
stressed (Ge:Ga) detectors w hereas the shorter wavelength 
channels were all unstressed ijClegg et al.|[l99^ . However, 
neither explanation has any corroborating evidence from 
other observations or from the literature. 
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Comparing the final throughput for the two FPs clearly 
shows that at short wavelengths it is advantageous in terms 
of signal-to-noise to use FPL non-prime data rather than 
FPS. The only cost is a slight reduction in spectral resolu- 
tion but this is outweighed by the signal-to-noise gain (see 
Sect. mi}. 



APPENDIX B: DATA 
APPENDIX C: OBSERVATIONS 
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Figure Bl. Pull LOS FP spectrum after co-adding prime and non-prime data using FPL for all LWS detectors (below 70 /im only 
non-prime data were used). The data are binned at 1/4 of the instrumental resolution element. Unassigned features detected at >3it are 
labelled 'U-line'. Two features axe labelled 'ARTIFACT' (at 57 ;um and 141.7 ^im). These are residual features due to the data reduction 
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Table CI. Log of all the Sgr B2 LOS observations and the corresponding LOl grating observation. The TDT (target dedicated time) 
number is an unique identifier for each observation. For the L03 data, only the wavelength range recorded on the prime detector is shown. 
In some observations the specified wavelength coverage was split between two prime detectors (this occurred at the edge of the detector 
nominal range). The length of the observation in seconds, the observation date, the number of repeated FP mini-scans and the LSR 
velocity correction applied are shown in the remaining columns. 



TDT 


Prime Wavelength 


FP 


Prime 


Length 


Date 


Repeated 
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Number 
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Detector 
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